ASYMPTOTICS OF SYMMETRIC POLYNOMIALS WITH 
APPLICATIONS TO STATISTICAL MECHANICS AND 
REPRESENTATION THEORY 



Abstract. We develop a new method for studying the asymptotics of 
symmetric polynomials of representation-theoretic origin as the num- 
ber of variables tends to infinity. Several applications of our method 
are presented: We prove a number of theorems concerning characters of 
infinite-dimensional unitary group and their g-deformations. We study 
the behavior of uniformly random lozenge tilings of large polygonal do- 
mains and find the GUE-eigenvalues distribution in the limit. We also 
investigate similar behavior for Alternating Sign Matrices ( equivalent ly, 
six-vertex model with domain wall boundary conditions). Finally, we 
compute the asymptotic expansion of certain observables in 0(n — 1) 
dense loop model. 
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1. Introduction 

1.1. Overview. In this article we study the asymptotic behavior of sym- 
metric functions of representation-theoretic origin, such as Schur rational 
functions or characters of symplectic or orthogonal groups, etc, as their num- 
ber of variables tends to infinity. In order to simplify the exposition we stick 
to Schur functions in the introduction where it is possible, but most of our 
results hold in a greater generality. 

The rational Schur function sx{xi, . . . ,Xn) is a symmetric Laurent poly- 
nomial in variables xi, Xn- They are parameterized by A^-tuples of 
integers A = (Ai>A2>--->AAr) (we call such A^-tuples signatures, they 
form the set GT^r) and are given by Weyl's character formula as 

r \ . , AT AlN 

det 

sxixi, . . .,xn) 



X-i+N-j 

X: 



Our aim is to study the asymptotic behavior of the normalized symmetric 
polynomials 

s\{xi, . . 



(1.1) Sx{xi,...,Xk]N,l) 



and also 



N 




„x Q, ^ sx{xi,...,Xk,l,q,q^,...,q^ ^ ^) 
(1.2) Sx{xi,...,Xk]N,q) = ^^3- , 

^ \\^ 1 ■ ■ ■ 1 q ) 

for some q > 0. Here A = A(A^) is allowed to vary with N, k is any fixed 
number and complex numbers, which may or may not vary 

together with A'^, depending on the context. Note that there are explicit 
expressions (Weyl's dimension formulas) for the denominators in formulas 



(1.1) and l\1.2h, therefore, their asymptotic behavior is straightforward. 



The asymptotic analysis of expressions (1.1), (1.2) is important because 
of the various applications in representation theory, statistical mechanics 
and probability, including: 

• For any k and any fixed xi, . . . , x^., such that \xi\ = 1, the con- 
vergence of Sx{xi, . . . , Xk; N,l) (from ( |1.1[ )) to some limit and the 
identification of this limit can be put in representation-theoretic 
framework as the approximation of indecomposable characters of 
the infinite-dimensional unitary group U{oo) by normalized charac- 
ters of the unitary groups U (N), the latter problem was first studied 
by Kerov and Vershik |VK| . 



The convergence of Sx{xi, . . . ,Xk; N,q) (from (1.2)) for any k and 
any fixed xi, . . . ,Xk is similarly related to the quantization of char- 
acters of U{oo), see [G] . 
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The asymptotic behavior of (1.1) can be put in the context of Ran- 



dom Matrix Theory as the study of the Harish-Chandra-Itzykson- 
Zuber integral 



(1.3) / exp{Tiace{AUBU-^))dU, 

Ju(N) 



'UiN) 

where A is a fixed Hermitian matrix of finite rank and B = B{N) 
is an N X N matrix changing in a regular way as — )• oo. In this 
formulation the problem was thoroughly studied by Guionnet and 
MaMa IGM,. 



A normalized Schur function (1.1) can be interpreted as the expec- 



tation of a certain observable in the probabilistic model of uniformly 
random lozenge tilings of planar domains. The asymptotic analysis 



of ( 1.1 ) as — oo with Xi = exp(yj/vA^) and fixed y^s gives a 
way to prove the local convergence of random tiling to a distribution 
of random matrix origin — the so-called GUE-minors process. In- 
formal argument explaining that such convergence should hold was 
suggested earlier by Okounkov and Reshetikhin in |ORlj . 
When A is a staircase Young diagram with 2A'" rows of lengths 



iV - 1, iV - 1, iV - 2, iV - 2, . . . , 1, 1, 0, 0, ( |LT| ) gives the expectation 
of a certain observable for the uniformly random configurations of 
the six-vertex model with domain wall boundary conditions (equiva- 
lently, Alternating Sign Matrices) . Asymptotic behavior as — >• cx) 
with Xi = exp{yi/y/N) and fixed yi gives a way to study the local 
limit of the six-vertex model with domain wall boundary conditions 
near the boundary. 



• For the same staircase A the expression involving (1.1) with k = 
4 and Schur polynomials replaced by the characters of symplectic 
group gives the mean of the boundary-to-boundary current for the 
completely packed 0{n = 1) dense loop model, see [GNPj . The 
asymptotics (now with fixed Xi, not depending on N) gives the limit 
behavior of this current, significant for the understanding of this 
model. 

In the present article we develop a new unified approach to study the 



asymptotics of normalized Schur functions (1.1), (1.2) (and also for more 
general symmetric functions like symplectic characters and polynomials cor- 
responding to the root system BCn), which gives a way to answer all of the 
above limit questions. There are 3 main ingredients of our method: 

(1) We find simple contour integral representations for the normalized 



Schur polynomials (1.1), (1.2) with k = 1, i.e. for 



-X sa(x,1,...,1) sx{x,l,q,...,q^ 2) 

and also for more general symmetric functions of representation- 
theoretic origin. 

(2) We study the asymptotics of the above contour integrals using the 
steepest descent method. 
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(3) We find formulas expressing (1.1), (1.2) as k x k determinants of 



expressions involving (1.4), and combining the asymptotics of these 



formulas with asymptotics of (1.4) compute limits of (1.1), (1.2). 



In the rest of the introduction we provide a more detailed description of 



our results. In Section 1.2 we briefly explain our methods. In Sections |1.3| 



1.4, 1.5, 1.6 1.7 we describe the applications of our method in asymptotic 
representation theory, probability and statistical mechanics. Finally, in Sec- 



tion 1.8 we compare our approach for studying the asymptotics of symmetric 
functions with other known methods. 

1.2. Our method. The main ingredient of our approach to the asymptotic 
analysis of symmetric functions is the following integral formula, which is 

(Ai > Aa > 



proved in Theorem 3.8 Let A 



sx{xi, 



• • > Aat), and let xi, 

) -^k ) 1 ) • • • 5 1) 



Xk 



be complex numbers. Denote 

Sxixi,...,Xk;N,l) 

sa(1, • • • , 1) 

with N — k Is in the numerator and Is in the denominator. 



Theorem 1.1 (Theorem 3.8). For any complex number x other than and 
1 we have 



(1.5) 



Sx{x;N,l) 



{N-l)l 1 



{x - 1)^-1 27ri Jc n£i(^ - (A,, + N-i)) 
where the contour C encloses all the singularities of the integrand. 



dz, 



We also prove various generalizations of formula (1.5): one can replace 
Is by the geometric series 1, q, q^, ■ ■ ■ (Theorem 3.6), Schur functions can be 
replaced with characters of symplectic group (Theorems |3.15 and 3.18) or, 
more, generally, with multivariate Jacobi polynomials (Theorem 3.22). In 
all these cases a normalized symmetric function is expressed as a contour 
integral with integrand being the product of elementary factors. The only 
exception is the most general case of Jacobi polynomials, where we have to 
use certain hypergeometric series. 

Recently (and independently of the present work) a formula similar to 



(1.5) for the characters of orthogonal groups 0{n) was found in [HJJ in the 
study of the mixing time of certain random walk on 0(n). A close relative 
of our formula (1.5) can be also found in Section 3 of |CPZj . 

Using formula (1.5) we apply tools from complex analysis, mainly the 



method of steepest descent, to compute the limit behavior of these normal- 
ized symmetric functions. Our main asymptotic results along these lines are 



summarized in Propositions 4.1 4.2, 4.3 for real x and in Propositions 4.4 



and 4.5 for complex x. 

The next important step is the formula expressing S\{xi, 
terms of S\{xi; N,l) which is proved in Theorem |3.7[ 



,Xk;N,l) in 



Theorem 1.2 (Theorem 3.7). We have 



(1.6) Sxixi,...,Xk;N,l) 



X 7 



n 



\N-k 
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X det 



1 k 



.. , \ X{Sx{xf,N,l) 



1 



(iV-1)! 



where is the differential operator x ^ . 



Formula ( 1.6 ) can again be generalized: Is can be replaced with geometric 
series 1, g, g"^, . . . (Theorem 3.5), Schur functions can be replaced with char- 
acters of the symplectic group (Theorems 3.14, 3.17) or, more, generally, 
with multivariate Jacobi polynomials (Theorem 3.21). In principle, formu- 
las similar to (1.6) can be found in the literature, see e.g. PPl Proposit ion 
6.2], [K^. 



The advantage of formula (1.6) is its relatively simple form, but it is not 
straightforward that this formula is suitable for the N ^ oo limit. However, 



we are able to rewrite this formula in a different form (see Proposition 3.9), 
from which this limit transition is immediate. Combining the limit formula 
with the asymptotic results for S\{x;N, 1) we get the full asymptotics for 
Sx{xi, . . . , Xk] N,l). As a side remark, since we deal with analytic functions 
and convergence in our formulas is always (at least locally) uniform, the 
differentiation in formula (1.6) does not introduce any problems. 



1.3. Application: asymptotic representation theory. Let U{N) de- 
note the group of all x unitary matrices. Embed U{N) into U{N + 1) 
as a subgroup acting on the space spanned by first N coordinate vectors 
and fixing A^ -|- 1st vector, and form the infinite-dimensional unitary group 
U{oo) as an inductive limit 



U{oo) = U U{N). 



N=l 



Recall that a (normalized) character of a group G is a continuous function 
Xia), 9 ^ G satisfying: 

(1) X is constant on conjugacy classes, i.e. x(o^fl~^) = x{t>), 



,9k}, 



is Hermitian 



(2) X is positive definite, i.e. the matrix 

non-negative definite, for any {gi, 

(3) X(e) = 1. 

An extreme character is an extreme point of the convex set of all charac- 
ters. If G is a compact group, then its extreme characters are normalized 
matrix traces of irreducible representations. It is a known fact (see e.g. the 
classical book of Weyl [W]) that irreducible representations of the unitary 
group U{N) are parameterized by signatures, and the value of the trace of 
the representation parameterized by A on a unitary matrix with eigenval- 
ues Ml, . . . , UN is sx{ui, . . . ,un)- Using these facts and applying the result 
above to U{N) we conclude that the normalized characters of U{N) are the 
functions 

sxjui,.. .,un) 
sa(1,...,1) 

For "big" groups such as [/(oo) the situation is more delicate. The study 
of characters of this group was initiated by Voiculescu |Voj in 1976 in con- 
nection with finite factor representations of U{oo). Voiculescu gave a list 
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of extreme characters, later independently Boyer [Boj and Vershik-Kerov 
|VK) discovered that the classification theorem for the characters of U{oo) 
follows from the result of Edrei jEdj on the characterization of totally posi- 
tive Toeplitz matrices. Nowadays, several other proofs of Voiculescu-Edrei 
classification theorem is known, see [00| . |BQ| . |P2| . The theorem itself 
reads: 

Theorem 1.3. The extreme characters of U{oo) are parameterized by the 
points 00 of the infinite- dimensional domain 

n C M^°°+^ = X ]R°^ X ]R°^ X R°° X R X R, 

where Vt is the set of sextuples 

uj = (a+,a-,/3+,/3-;(5+,5-) 

such that 

a± = (af > > • • • > 0) € R°°, = > /Ss^ > • • • > 0) G R~, 

oo 

5^(«f +/3f)<<5^, /3+ + /3r<l. 

i=l 

The corresponding extreme character is given by the formula 
(1-7) 

X'^^\U)= TT ^5+(u-l)+5-(u-^-l) TT 1 + f^tiu - 1) 1 + PAu-^ - 1) 

.espe!tL(c/) li 1 - ""ti^ - 1) 1 - - 1) 

where 

oo 
i=l 

Our interest in characters is based on the following fact. 

Proposition 1.4. Every extreme normalized character x ofU{(x>) is a uni- 
form limit of extreme characters ofU{N). In other words, for every x there 
exists a sequence \{N) G GTat such that for every k 

x(ni,...,'Ufc,l,...) = lim Sx{ui, . . . ,Uk; N,l) 

uniformly on the torus (Si)^, where Si = {n G C : |u| = 1}. 

In the context of representation theory of C/(oo) this statement was first 
observed by Kerov and Vershik [VK| . However, this is just a particular case 
of a very general convex analysis theorem which was reproved many times 
in various contexts (see e.g. [Vj, [OOj . |DF2j ). 

The above proposition raises the question which sequences of characters 
of U{N) approximate characters of C/(oo). Solution to this problem was 
given by Kerov and Vershik |VK) . 

Let // be a Young diagram with row lengths fii, column lengths fi[ and 
whose length of main diagonal is d. Introduce modified Frobenius coordi- 
nates: 

Pi = Hi- i + 1/2, qi = fi'i-i + 1/2, i = l,...,d. 
Note that ^^^=1 Pi + Qi = 1/^1 • 
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Given a signature A G GT^v, we associate two Young diagrams and 
A~ to it: The row lengths of A^ are the positive Aj's, while the row lengths 
of A~ are minus the negative ones. In this way we get two sets of modified 
Frobenius coordinates: pf ,qf , i = I, . . . ,d~^ and p~ ,q~ , i = 1, . . . ,d~ . 

Theorem 1.5 ([VK], [00], [HO], [El]). Let uj = (a±,/3±;(5±) and suppose 
that the sequence \{N) G GTjy is such that 

ptiN)/N^at, pr{N)/N^ar, q+(N)/N^Pt, q;{N)/N^/3t, 
|A+|/7V^(5+, \X-\/N^6-. 

Then for every k 

x'^(ni,...,Ufc,l,...) = lim Sxtmiui, ■ ■ ■ ,Uk; N,l) 

uniformly on the torus (Si)^. 

Theorem |1.5| is an immediate corollary of our results on asymptotics of 



normalized Schur polynomials, and a new short proof is given in Section 6.1 
Note the remarkable multiplicativity of Voiculescu-Edrei formula for the 
characters of U{oo): the value of a character on a given matrix (element 
of U{oo)) is expressed as a product of the values of a single function at 
each of its eigenvalues. There exists an independent representation-theoretic 
argument explaining this multiplicativity. Clearly, no such multiplicativity 
exists for finite N, i.e. for the characters of U{N). However, we claim 
that the formula ( |1.6[ ) should be viewed as a manifestation of approximate 
multiplicativity for (normalized) characters of U{N). To explain this point 
of view we start from A; = 2, in this case (|1.6[) simplifies to 



Sx{x,y;N,l) = Sx{x;N,l)Sx{y;N,l) 

- l)(y - 1) (xl - yly)[Sx{x; N, l)Sx{y; N, 1)] 



+ 



— 1 y — X 



More generally Proposition 3.9 claims that for any k formula (1.6) implies 
that, informally, 

Sa(xi, . . . , Xfc; iV, 1) = Sx{xi;N, 1) • • • 5a (x^; N, 1) + 0(l/iV), 



therefore, ( 1.6 ) states that normalized characters of U (N) are approximately 
multiplicative and they become multiplicative as — )• cxd. This is somehow 
similar to the work of Diaconis and Freedman [DFJ on finite exchangeable 
sequences. In particular, in the same way as results of [DF\ immediately 
imply de Finetti's theorem (see e.g. [A]), our results immediately imply the 
multiplicativity of characters of U{oo). 

In [Gj a g-deformation of the notion of character of U{oo) was suggested. 



Analogously to Proposition 1.4, a g-character is a limit of Schur functions, 
but with different normalization. This time the sequence A(A^) should be 
such that for every k 

, . Sx(N){xi,---,Xk,q '',q ,---,q' 

^ ' SxiN){l,q-\...,q'-'') 

converges uniformly on the set {(xi, . . . ,Xk) G C'^ | |xi| = q^~^}. An ana- 
logue of Theorem |1.5| is the following one: 
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Theorem 1.6 (^). Let < q < 1. Extreme q-characters of U{oo) are 
parameterized by the points of set M of all non- decreasing sequences of in- 
tegers: 

AA = {z^i < 1/2 < 1^3 < •••} C 
Suppose that a sequence X{N) G GTjy is such that for any j > 
(1.9) lim Xn+i-j{N) 



then for every k 



(1.10) 



-k-l 



converges uniformly on the set {(xi, 
limits define the q-character of U (oo 



5A{JV)(l,'?-^•••,'?l-^ 



q *} and these 



Using the g-analogues of formulas (1.5) and (1.6) we give in Section 6.2 



a short proof of the second part of Theorem 1.6 see Theorem 6.5 This 
should be compared with [G], where the proof of the same statement was 
quite involved. We go beyond the results of [G], give new formulas for 
the g-characters and explain what property replaces the multiplicativity of 



Voiculescu-Edrei characters given in Theorem 1.3 



1.4. Application: random lozenge tilings. Consider a tiling of a domain 
drawn on the regular triangular lattice of the kind shown at Figure [l]with 
rhombi of 3 types, where each rhombus is a union of 2 elementary triangles. 
Such rhombi are usually called lozenges and they are shown at Figure [2] The 
configuration of the domain is encoded by the number N which is its width 
and N integers /ii > ^2 > • • • > tJ-N which are the positions of horizontal 
lozenges sticking out of the right boundary. If we write //j = Xi + N — i, then 
A is a signature of size A^, see left panel of Figure [!} Due to combinatorial 
constraints the tilings of such domain are in correspondence with tilings of 
a certain polygonal domain, as shown on the right panel of Figure [T} Let 
Q,x denote the domain encoded by a signature A. 

It is well-known that each lozenge tiling can be identified with a stepped 
surface in (the three types of lozenges correspond to the three slopes 
of this surface) and with a perfect matching of a subgraph of a hexagonal 
lattice, see e.g. |Kej . Note that there are finitely many tilings of Q\ and 
let T;^ denote a uniformly random lozenge tiling of 0,x. The interest in 
lozenge tilings is caused by their remarkable asymptotic behavior. When 
is large the rescaled stepped surface corresponding to T;s^ concentrates near a 
deterministic limit shape. In fact, this is true also for more general domains, 
see |CKP| . One feature of the limit shape is the formation of so-called 
frozen regions; in terms of tilings, these are the regions where asymptotically 
with high probability only single type of lozenges is observed. This effect is 
visualized in Figure [3j where a sample from the uniform measure on tilings of 
the simplest tilable domain — hexagon — is shown. It is known that in this 
case the boundary of the frozen region is the inscribed ellipse, see |CLPj . 
for more general polygonal domains the frozen boundary is an inscribed 
algebraic curve, see |K0] and also jPlj . 
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Figure 1. Lozenge tiling of the domain encoded by signa- 
ture A (left panel) and of corresponding polygonal domain 
(right panel). 



Figure 2. The 3 types of lozenges, the middle one is called 
"horizontal" . 

In this article we study the local behavior of lozenge tiling near a turn- 
ing point of the frozen boundary, which is the point where the boundary of 
the frozen region touches (and is tangent to) the boundary of the domain. 
Okounkov and Reshetikhin gave in |ORlj a non-rigorous argument explain- 
ing that the scaling limit of a tiling in such situation should be governed 
by the so-called GUE-minors process (introduced and studied by Barysh- 
nikov [Bar] and Johansson-Nordenstam [JNj), which is the joint distribution 
of the eigenvalues of a Gaussian Unitary Ensemble (GUE-)random matrix 
(i.e. Hermitian matrix with independent Gaussian entries) and of its top- 
left corner square submatrices. In one model of tilings of infinite polygonal 
domains, the proof of the convergence can be based on the determinantal 
structure of the correlation functions of the model and on the double-integral 
representation for the correlation kernel and it was given in [ORT] . Another 
rigorous argument, related to the asymptotics of orthogonal polynomials ex- 
ists for the lozenge tilings of hexagon (as in Figure |3]), see |JN) . [N] . 

Given T\ let I'l > 1^2 > ■ ■ ■ > I'k be the horizontal lozenges at the kth 
vertical line from the left. (Horizontal lozenges are shown in blue in the left 



10 



VADIM GORIN AND GRETA PANOVA 




Figure 3. A sample from uniform distribution on tilings of 
40 X 50 X 50 hexagon and corresponding theoretical frozen 
boundary. The three types of lozenges are shown in three 
distinct colors. 



panel of Figure [T]) We set Ui = Ki -\- k — i and denote the resulting random 
signature k of size k as T^. Further, let GUEfc denote the distribution 
of k (ordered) eigenvalues of a random Hermitian matrix from a Gaussian 
Unitary Ensemble. 



Theorem 1.7 (Theorem [5T|). Let X{N) G GTn, N = 1,2,... be a se- 
quence of signatures. Suppose that there exist a non-constant piecewise- 
differentiable weakly decreasing function f{t) such that 



E 

i=l 



urn 



N 



as N ^ oo and also supj ^ |Ai(A^)/A^| < oo. Then for every k as N ^ oo 
we have 

T^,^. - NE{f) 
^/NSU) 

in the sense of weak convergence, where 

E{f) = C f{t)dt, S{f) = C f{tfdt - E{ff + /' /(t)(l - 2t)dt. 
Jo Jo Jo 

Corollary 1.8 (Corollary|5.2[). Under the same assumptions as in Theorem 



1.1 the (resettled) joint distribution of k(k + l)/2 horizontal lozenges on the 



left k lines weakly converges to the joint distribution of the eigenvalues of 
the k top-left corners of a k x k matrix from a GUE ensemble. 

Note that, in principle, our domains may approximate a non-polygonal 
limit domain as — )• oo, thus, the results of |KOj describing the limit 
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shape in terms of algebraic curves are not applicable here and not much is 
known about the exact shape of the frozen boundary. In particular, even the 
explicit expression for the coordinate of the point where the frozen boundary 



touches the left boundary (which we get as a side result of Theorem 1.7) 
seems not to be present in the literature. 



Our approach to the proof of Theorem 1.7 is the following: We express 
the expectations of certain observables of uniformly random lozenge tilings 
through normalized Schur polynomials S\ and investigate the asymptotics 
of these polynomials. In this case we prove and use the following asymptotic 



expansion (given in Proposition 4.3 and Proposition 5.8) 

Sx{e^,...,e^;N,l) 

= exp (^VNE{f){hi + ... + hk) + ls{f){hj + ... + hl) + o(i; 

We believe that our approach can be extended to a natural q'-deformation 
of uniform measure, which assigns the weight q"""^ to lozenge tiling with 
volume vol below the corresponding stepped surface; and also to lozenge 
tilings with axial symmetry, as in |FNj . |BKj . In the latter cases the Schur 
polynomials are replaced with characters of orthogonal or symplectic groups 
and the limit object also changes. We postpone the thorough study of these 
cases to a future publication. 

We note that there might be another approach to the proof of Theorem 



1.7 Recently there was progress in understanding random tilings of polygo- 
nal domains, Petrov found double integral representations for the correlation 
kernel describing the local structure of tilings of a wide class of polygonal 
domains, see |Plj (and also |Me] for a similar result in context of random 
matrices). Starting from these formulas, one could try to prove the GUE- 
minors asymptotics along the lines of [ORlj . 

1.5. Application: six vertex model and random ASMs. An Alter- 
nating Sign Matrix of size is a x matrix whose entries are either 
0, 1 or —1, such that the sum along every row and column is 1 and, more- 
over, along each row and each column the nonzero entries alternate in sign. 
Alternating Sign Matrices are in bijection with configurations of the six- 
vertex model with domain- wall boundary conditions as shown at Figure l4l 



more details on this bijection are given in Section 5.2 A good review of the 



six-vertex model can be found e.g. in the book [Baxj by Baxter. 

Interest in ASMs from combinatorial perspective emerged since their dis- 
covery in connection with Dodgson condensation algorithm for determinant 
evaluations. Initially, questions concerned enumeration problems, for in- 
stance, finding the total number of ASMs of given size n (this was the long- 
standing ASM conjecture proved by Zeilberger fZe^ and Kuperberg [Kuj . 
the full story can be found in the Bressoud's book [Br]). Physicists' inter- 
est stems from the fact that ASMs are in one-to-one bijection with con- 
figurations of the six-vertex model. Many questions on ASMs still remain 
open. Examples of recent breakthroughs include the Razumov-Stroganov 
|RSj conjecture relating ASMs to yet another model of statistical mechan- 
ics (so-called 0(1) loop model), which was finally proved very recently by 
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1 
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1 0\ 

-1 1 





1 oj 




Figure 4. An alternating sign matrix of size 5 and the cor- 
responding configuration of 6-vertex model with domain wall 
boundary condition. 

Cantinni and Spontiello [CS], and the still open question on a bijective proof 
of the fact that Totally Symmetric Self-Complementary Plane Partitions and 
ASMs are equinumerous. A brief up-to-date introduction to the subject can 
be found e.g. in |BFZj . 

Our interest in ASMs and the six-vertex model is probabilistic. We would 
like to know how a uniformly random ASM of size n looks like when n is 
large. Conjecturally, the features of this model should be similar to those 
of lozenge tilings: we expect the formation of a limit shape and various 
connections with random matrices. The existence and properties of the limit 
shape were studied by Colomo and Pronko [CP3j . however their arguments 
are physical, while mathematical proof is yet unavailable. 

In the present article we prove a partial result toward the following con- 
jecture. 

Conjecture 1.9. Fix any k. As n ^ oo the probability that the number of 
(— l)s in the first k rows of a uniformly random ASM of size n is maximal 
(i.e. there is one ( — 1) in second row, two (— l)s in third row, etc) tends to 
1, and, thus, Is in first k rows are interlacing. After proper centering and 
resettling, the distribution of the positions of Is tends to the GUE-minors 
process as n ^ oo. 

Let ^ki'iT') denote the sum of coordinates of Is minus the sum of coordi- 
nates of ( — l)s in the kth row of the uniformly random ASM of size n. We 
prove that the centered and rescaled random variables ^fc(n) converge to 
the collection of i.i.d. Gaussian random variables as n — t- oo. 



Theorem 1.10 (Theorem 5.9). For any fixed k the random variable 

^fe(n)-n/2 



weakly converges to the normal random variable N{0, yS/S). 
Moreover, the joint distribution of any collection of such variables converges 
to the distribution of independent normal random variables N{0, y^3/8). 
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Remark. We also prove a bit stronger statement, see Theorem 5.9 for the 
details. 

Indeed, if the latter 



Note that Theorem 1.10 agrees with Conjecture 1.9 



holds, then ^^kin) converges to the difference of the sums of the eigenvalues 
oi akxk GUE-random matrix and of its {k — l)x(k — 1) top left submatrix. 
But these sums are the same as the traces of the corresponding matrices, 
therefore, the difference of sums equals the bottom right matrix element of 
the k X k matrix, which is a Gaussian random variable by the definition of 
GUE. 

Our proof of Theorem 1.10 has two components. First, a result of Okada 
|Qkj . based on earlier work of Izergin and Korepin (ij, |Kor] . shows that 
sums of certain quantities over all ASMs can be expressed through Schur 
polynomials (in an equivalent form this was also shown by Stroganov [St]). 
Second, our method gives the asymptotic analysis of these polynomials. 
In fact, we claim that Theorem 1.10 together with additional probabilistic 



argument implies Conjecture 1.9 However, this argument is unrelated to 
the asymptotics of symmetric polynomials and, thus, is left out of the scope 
of the (already long) present paper and is postponed to a future publication. 

In the literature one can find another probability measure on ASMs as- 
signing the weight 2"^ to the matrix with m ones. For this measure there 
are many rigorous mathematical results, due to the connection to the uni- 
form measure on domino tilings of the Aztec diamond, see [EKLPj . |FS| . 
The latter measure can be viewed as a determinantal point process, which 
gives tools for its analysis. An analogue of Conjecture 1.9 for the tilings of 
Aztec diamond was proved by Johansson and Nordenstam |JN) . 

In regard to the combinatorial questions on ASMs, we note that there 
has been interest in refined enumerations of Alternating Sign Matrices, i.e. 
counting the number of ASMs with fixed positions of Is along the boundary. 
In particular, Colomo-Pronko |CPlj . |CP2j . Behrend |Bej and Ayyer-Romik 
|AR ] found formulas relating A;~refined enumerations to 1-refined enumera- 
tions for ASMs. Some of these formulas are closely related to particular cases 
of our multivariate formulas (Theorem 3.7) for staircase Young diagrams. 



1.6. Application: 0{n = 1)— loop model. Recently found parafermionic 
observables in the so-called completely packed 0{n = 1) dense loop model 
in a strip are also simply related to symmetric polynomials, see [GNPj . The 
0{n = 1) dense loop model is one of the representations of the percola- 
tion model on the square lattice. For the critical percolation models similar 
observables and their asymptotic behavior were studied (see e.g. |iSmj ) . how- 
ever, the methods involved are usually completely different from ours. 

A configuration of the 0{n = 1) loop model in a vertical strip consists 
of two parts: a tiling of the strip on a square grid of width L and infinite 
height with squares of two types shown in Figure [5] (left panel) , and a choice 
of one of the two types of boundary conditions for each 1x2 segment along 
each of the vertical boundaries of the strip; the types appearing at the left 
boundary are shown in Figure [5] (right panel) . Let ~I l denote the set of 
all configurations of the model in the strip of width L. An element of "Ig 
is shown in Figure [5| Note that the arcs drawn on squares and boundary 
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segments form closed loops and paths joining the boundaries. Therefore, 
the elements of "l^, have an interpretation as collections of non- intersecting 
paths and closed loops. 





Figure 5. Left panel: the two types of squares. Right 
panel: the two types of boundary conditions. 
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Figure 6. A particular configuration of the dense loop 
model showing a path passing between two vertically adja- 
cent points X and y. 



In the simplest homogeneous case a probability distribution on "1^ is de- 
fined by declaring the choice of one of the two types of squares to be an 
independent Bernoulli random variable for each square of the strip and for 
each segment of the boundary. I.e. for each square of the strip we flip an 
unbiased coin to choose one of the two types of squares (shown in Figure 
[5]) and similarly for the boundary conditions. More generally, the type of a 
square is chosen using a (possibly signed or even complex) weight defined as 
a certain function of its horizontal coordinate and depending on L parame- 
ters zi, . . . ,zi; two other parameters Cij C2 control the probabilities of the 
boundary conditions and, using a parameter q, the whole configuration is 
further weighted by its number of closed loops. We refer the reader to [GNP] 
and references therein for the exact dependence of weights on the parameters 
of the model and for the explanation of the choices of parameters. 

Fix two points x and y and consider a configuration w G "Ix,. There are 
finitely many paths passing between x and y. For each such path r we 
define the current c(r) as if r is a closed loop or joins points of the same 
boundary; 1 if r joins the two boundaries and x lies above r; — 1 if r joins 
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the two boundaries and x hes below r. The total current C^'^(a;) is the sum 
of c(t) over all paths passing between x and y. The mean total current F^'^ 
is defined as the expectation of C^'^. 

Two important properties of F^'^ are skew-symmetry 

px,y _ _py,x 

and additivity 

pXi,Xz _ pXi,X2 _|_ pX2,X3 

These properties allow to express as a sum of several instances of the 

mean total current between two horizontally adjacent points 

_p(*.i):(*.i+i) 

and the mean total current between two vertically adjacent points 

_p(i.«),{i+i>«)_ 

The authors of [GNP] present a formula for and F^^'^^'^^+i'*) 

which, based on certain assumptions, expresses them through the symplec- 
tic characters Xx^i^l 4X1,(1) where = ( [^J , [^J , • ■ • , 1, 0, 0). 
The precise relationship is given in Section |5.3[ Our approach allows us 
to compute the asymptotic behavior of the formulas of |GNP] as the lat- 



tice width L — )• cx), see Theorem 5.12 In particular, we prove that the 



leading term in the asymptotic expansion is independent of the boundary 
parameters Ci, C2- 

This problem was presented to the authors by de Gier [Gij. |GPj during 
the program "Random Spatial Processes" at MSRI, Berkeley. 

1.7. Application: matrix integrals. Let A and B be two N x N Hermit- 
ian matrices with eigenvalues ai, . . . , oat and bi, . . . ,bi^, respectively. The 
Harish-Chandra formula |H1| . |H2j (sometimes known also as Itzykson- 
Zuber [IZ] formula in physics literature) is the following evaluation of the 
integral over the unitary group: 
(1.11) 

r , detij=i,,..,Ar(exp(ai6j) ) 

Ju{N) Ui<j (fli - aj ) Ui<j (bi - bj ) f^. 

where the integration is with respect to the normalized Haar measure on 
the unitary group U{N). Comparing ( |1.11[ ) with the definition of Schur 
polynomials and using Weyl's dimension formula 

.>(i.....i)^n "-"':f'~^' . 

i<j 

we observe that when bj = Xj + N — j the above matrix integral is the 
normalized Schur polynomial times explicit product, i.e. 

SA(e'^S...,e'^") 



— TT 

sa(1,.--,1) fj^. ai-aj 

Guionnet and Maida studied (after some previous results in the physics 
literature, see |GMj and references therein) the asymptotics of the above 
integral as — )• oo when the rank of A is finite and does not depend on 
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A^. This is precisely the asymptotics of (1.1). Therefore, our methods (in 



particular, Propositions 4.1 , 4.2 4.4 ) give a new proof of some of the results 
of [GM] . In the context of random matrices the asymptotics of this integral 
in the case when rank of A grows as the size of A grows was also studied, 
see e.g. |GZ] . |CSj . However, currently we are unable to use our methods 
for this case. 

1.8. Comparison with other approaches. Since asymptotics of symmet- 
ric polynomials as the number of variables tends to infinity already appeared 
in various contexts in the literature, it makes sense to compare our approach 
to the ones used before. 

In the context of asymptotic representation theory the known approach 
(see |VKj . |Q0] . |002] . [G]) is to use the so-called binomial formulas. In 
the simplest case of Schur polynomials such formula reads as 

(1.12) Sxil + xi,...,l + Xk;N,l) = Y^ s^ixi, XkHp, A, N), 

where the sum is taken over all Young diagrams ^ with at most k rows, 
and c(/i. A, A^) are certain (explicit) coefficients. In the asymptotic regime of 



Theorem |1.5| the convergence of the left side of ( 1.12 ) implies the convergence 
of numbers c(/i. A, N) to finite limits as — )• oo. Studying the possible 
asymptotic behavior of these numbers one proves the limit theorems for 
normalized Schur polynomials. 

Another approach uses the decomposition 

(1.13) 5'a(2;i, . . . ,Xk]N,l) = ^ S^{xi, ...,Xk;k, l)d{fi, A, A^), 

where the sum is taken over all signatures of length k. Recently in |BOj and 
|P2j k X k determinantal formulas were found for the coefficients d{n, A, A'^). 
Again, these formulas allow the asymptotic analysis which leads to the limit 
theorems for normalized Schur polynomials. 



The asymptotic regime of Theorem 1.5 is distinguished by the fact that 
|Ai(A^)|/A^ is bounded as A^ — )■ oo. This no longer holds when one studies 
asymptotics of lozenge tilings, ASMs, or 0{n = 1) loop model. As far as 



authors know, in the latter limit regime neither formulas (1.12) nor (1.13) 
give simple ways to compute the asymptotics. The reason for that is the 
fact that for any fixed jj, both c(/x. A, N) and (i(/x. A, N) would converge to 
zero as A^ — )• cxD and more delicate analysis would be required to reconstruct 
the asymptotics of normalized Schur polynomials. 



Yet another, but similar approach to the proof of Theorem 1.5 was used 
in |Bo2j but, as far as authors know, it also does not extend to the regime 
we need for other applications. 

On the other hand the random-matrix asymptotic regime of [GM] is sim- 
ilar to the one we need for studying lozenge tilings, ASMs, or 0{n = 1) 
loop model. The approach of |GM| is based on the matrix model and 
the proofs rely on large deviations for Gaussian random variables. How- 
ever, it seems that the results of |GMj do not suffice to obtain our appli- 
cations: for k > 1 only the first order asymptotics (which is the limit of 
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lji{S\{xi, . . . , Xk] N, 1))/N) was obtained in |GMj . while our applications re- 
quire more delicate analysis. It also seems that the results of [GM] ( even for 
A; = 1) cannot be applied in the framework of the representation theoretic 



regime of Theorem 1.5 
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2. Definitions and problem setup 

In this section we set up notations and introduce the symmetric functions 
of our interest. 

A partition (or a Young diagram) A is a collection of non-negative numbers 
Ai > A2 > . . . , such that Aj < 00. The numbers Aj are row lengths of A, 
and the numbers A^ = \{j : Xj > i}\ are column lengths of A. 

More generally a signature A of size is an A^-tuple of integers Ai > 
A2 > • • • > Aat. The set of all signatures of size is denoted GTjy. It is 
also convenient to introduce strict signatures, which are A^-tuples satisfying 
strict inequalities Ai > A2 > • • • > Xjy; they from the set GT^r. We are going 
to use the following identification between elements of GTjv and GT^v: 

GTn3X< — > X + 6n = I^^G^Tn, fii = Xi + N-i, 
where we set 6n = {N — 1, N — 2, . . . ,1). The subset of GT^r (GT^v) of all 



signatures (strict signatures) with non-negative coordinates is denoted GTf 



(GT^). 

One of the main objects of study in this paper are the rational Schur func- 
tions, which originate as the characters of the irreducible representations of 
the unitary group U{N) ( equivalent ly, of irreducible rational representa- 
tions of the general linear group GL{N)). Irreducible representations are 
parameterized by elements of GT^r, which are identified with the dominant 
weights, see e.g. jWj or |Zhj . The value of the character of the irreducible 
representation V\ indexed by A G GT^v, on a unitary matrix with eigenvalues 
til, . . . , Mat is given by the Schur function, 

N 

det 

(2.1) sx{ui,...,un) 



X,+N~j 



U^<jiui-Uj) 

which is a symmetric Laurent polynomial in ui, . . . ,un- The denominator 



in (2.1) is the Vandermonde determinant and we denote it through A: 
A(ui, at) = det lij V . = \\{ui-Uj). 
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When the numbers Ui form a geometric progression, the determinant in (2.1 ) 
can be evaluated expUcitly as 

„\,+N~i _ \j+N-j 

(2.2) ■^"-'j^n ' ■ 

i<j 

In particular, sending g — )• 1 we get 

(2.3) .,(1^)= n ^^'-'\-^^^-'\ 

l<i<j<N •' 

where we used the notation 

1^ = 




The identity (2.3) gives the dimension of Vx and is known as the Weyl's 
dimension formula. 

In what follows we intensively use the normalized versions of Schur func- 
tions: 

Sx{xi,...,Xk;N,q) = — -^^3- , 

in particular, 

Sx{xi,...,Xk,l^~'') 



Sxixi, ...,Xk;N, 1) 



sa(I^) 

The Schur functions are characters of type A (according to the classifi- 
cation of root systems), their analogues for other types are related to the 
multivariate Jacobi polynomials. 

For a,b > —1 and m = 0, 1, 2 . . . let pmix; a, b) denote the classical Jacobi 
polynomials orthogonal with respect to the weig ht (1 - x)"(l + xf on the 
interval [—1,1], see e.g. |Erl] . |KoSw| . We use the normalization of |Erlj . 
thus, the polynomials can be related to the Gauss hypergeometric function 
2-^1: 

fffi -\- Q, -\- \^ I 1 — X 

pmix; a, b) = — — T:p7 — — V 2F1 -m, m + a + b + l,a + I; — — 

i{m+ Ijl (a -|- 1) \ 2 

For any strict signature A E GT^y set 

m / n det[pA.(xj;a,6)]f.^i 

^ix{xi, . . . ,XN;a,b) = — . 

A(xi, . . .,xn) 

and for any (non-strict) A € GT^ define 

(2.4) ...,z^;a,b)= cx^x+S • • • , «, , 

where cx is a constant chosen so that the leading coefficient of ^x is 1- 
The polynomials ^x are (a particular case of) BCn multivariate Jacobi 
polynomials, see e.g. |0Q2j and also |HS] . |M2j . |Kooj . We also use their 
normalized versions 

(2.5) Jx{zi,. . . ,Zk;N,a,b) = . 



ASYMPTOTICS OF SYMMETRIC POLYNOMIALS 



19 



Again, there is an explicit formula for the denominator in (2.5) and also for 
its g-version. For special values of parameters a and b the functions Jx can be 
identified with spherical functions of classical Riemannian symmetric spaces 
of compact type, in particular, with normalized characters of orthogonal and 
symplectic groups, see e.g. |Q02l Section 6]. 

Let us give more details on the latter case of the symplectic group Sp{2N), 
as we need it for one of our applications. This case corresponds to a = 6 = 
1/2 and here the formulas can be simplified. 

The value of character of irreducible representation of Sp{2N) parameter- 
ized by A G GT]^ on symplectic matrix with eigenvalues xi, x^^ , . . . , xn,x^^ 
is given by (see e.g. jWj, |Zhj ) 

TV 



det 



Xx{xi, . ■ ■,xn) 



\j+N+l-j 



-\j-N+l-j 



det 



N+l-j -Af-l+j 



N 



The denominator in the last formula can be evaluated explicitly and we 
denote it Ag 



(2.6) As(xi, . . . ,XAr) = det 



N-j+l 



N 



llixi-x-^) ll{x^+xr^-{xj+x-^)) 



UiKji^i - Xj){XiXj - 1) ll-{xf - 1) 



i<j 



(xi ■ • -Xn)" 



The normalized symplectic character is then defined as 



Xa(xi, ...,Xk;N, q) 



in particular 



Xxixi,...,Xk;N,l) 



Xx{xi, . . . ,Xfc,g, . . . , 
X\{xi,...,Xk,l^~'') 



Xa(I^) 

and both denominators again admit explicit formulas. 

In most general terms, in the present article we study the symmetric 
functions S\, Jx, Xx, their asymptotics as N ^ oo and its applications. 



Some further notations. 

We intensively use the (^-algebra notations 



1 



ml 



1 ' 



ni 

m=l 



and g-Pochhammer symbol 



fc-i 



i=0 



Since there are lots of summations and products in the text where i plays 
the role of the index, we write i for the imaginary unit to avoid the confusion. 
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3. Integral and multivariate formulas 



In this section we derive integral formulas for normalized characters of 
one variable and also express the multivariate normalized characters as de- 
terminants of differential (or, sometimes, difference) operators applied to 
the product of the single variable normalized characters. 

We first exhibit some general formulas, which we later specialize to the 
cases of Schur functions, symplectic characters and multivariate Jacobi poly- 
nomials. 

3.1. General approach. 

Definition 3.1. For a given sequence of numbers 6 = (61,62, ■■ ■), a col- 
lection of functions . . . , xat)}, N = 1,2, . . . , fj, G GTn (orGTj^) 
is called a class of determinantal symmetric functions with parameter 6, if 
there exist functions a{u), f3{u), g{u,v), numbers cn, and linear operator T 
such that for all N and fi we have 

(1) 



A^{xi, . . .,xn) 



det[g{xi; 
A(xi, . . .,xn) 



(2) 



N 



Af,{6i, . . . ,6n) = cnY\_ - ainj)), 

i=l i<j 

(3) g{x; m) (m G Z for the case of GT and m G Z>o for the case GT ) 
are eigenf unctions ofT acting on x with eigenvalues a{m), i.e. 

T{g{x,m)) = a{m)g{x,m). 

(4) a'{m) 7^ for all m as above. 



Proposition 3.2. For A^{xi, . . . ,xn) , as in Definition 3.1 we have the 
following formula 



(3.1) 



Ap,{xi, . . . ,Xk,6i, . . . ,di^_k) CN-k 



det 



A^{6i, . . . , 6n) 

k 



CN 



k N-k 

nn 

i=i j=i 



1 



-1 



A(xi, ...,Xk) 



n 

i=l 



A^jxj, 61, ... , 6n-i) 
Afj,{6i, . . . , 6n) 



N-l 



^ CAT 
CN-l 



where Ti is operator T acting on variable re,. 

Remark. Since operators Tj commute, we have 



det 



k 



i<j 



We also note that some of the denominators in (3.1) can be grouped in the 
compact form 

k N-k 

n n (^^i - %)A(xi, ...,xk) 
i=i j=i 



A{xi, . . . ,Xk,6i, . . . 



, ON-kj 



A(0i, . . . , Gj^^k) 
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Moreover, in our applications, the coefficients Cm will be inversely propor- 
tional to A{Oi, . . . ,9m), so we will be able to write alternative formulas 
where the prefactors are simple ratios of Vandermondes. 



Proof of Proposition 3.2. We will compute the determinant from property 
(1) of ^ by summing over all A; x A; minors in the first k columns. The rows 
will be indexed by / = {ii < ^2 < • • • < ik} and /x/ = (/Uj^ , . . . , /^jj.)- I'^ will 
be the complement of I in [l..n]. We have 



(3.2) 



'1, . . . ,UN~k} 



ntinf=?(^. 



•A.^{Oi, ■ ■ ■ , On) 

E {-l)^'^'^'-'^A,,ix,,...,Xk) 

I={ii<i2<--<ik} 

For each set / we have 



A 



flic 1 



9N-k) 



Afj,{6i, . . . , 6n) 



(3.3) 



■4^,,(gi,...,gjv-fc) CN ^ Uiei- ^(.f^i) Ui<j;i,j(^Aa(.f^i) - 

Af,{Ou...,eN) CN-k llZiPif^i)Ui<i<j<Ni»if^i) -(^(.f^j)) 



n^nn 



n 



^eI t^ti - "(^^-^ ^0,f^,^<J "^^^^^ " "^^^'^ 

-[-[- 1 -[-[- -[-[- 1 UieiUi<j,jeiHf^i)-(^i^j)) 



■1) 



i-1 



i{a{ni) - a{nj)) 



We also have that 



(3-4) Jl {a{fii) -a{fij))A^j{xi,...,Xk)A{xi,...,Xk) 

-ik k 



det 



a(/i. 



ley 



^=1 



E(-irdet 



= Y.{-l)^det 



det 



/-I 



Combining (3.2), (3.3) and (3.4) we get 

k N-k 



(3.5) 



A^{xi, . . . ,Xk,Oi, . . . , 9N-k) 



A^{9i, . . . , 9n) 



det [ly 



1 i,j=l 



A(xi, . . . ,a;fc) ^_ 



nn 

i=i j=i 

E En 

■■{ii<i2<-<ik} (TGS{k) i 



^ CN 

CN-k 
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Note that double summation in the last formula is a summation over all (or- 
dered) collections of distinct numbers. We can also include into the sum the 
terms where some indices ii coincide, since application of the Vandermonde 
of linear operators annihilates such terms. Therefore, (|3.2|) equals 



A(xi, ...,Xk) f}^^^ - a{^ij)) 



N 



When k = 1 the operators and the product over i disappear, 
so we see that the remaining sum is exactly the univariate ratio 

A^{xe,Oi, . . . ,9n-i) ^TT , CAT 



-^tiiOi, . . . , On) 
mula. 



X 7V-1 



CN-l 



and we obtain the desired for- 



□ 



Proposition 3.3. Under the assumptions of Definition \3.1\ we have the 
following integral formula for the normalized univariate Afj, 



(3.6) 



6n- 



CN-l 
CN 



N-1 

n 

i=l 



1 



1 



Oi / 27ri 



g{x;z)a'{z) 



c /3(2)n^Ii(«(^) - "(^i)) 



-dz. 



Here the contour C includes only the poles of the integrand at z = fn, i 
1 A^. 



Proof. As a byproduct in the proof of Proposition |3.2| we obtained the fol- 
lowing formula: 



(3.7) 



Af,{x,9i, . . .,9n-i) 
Afj,{9i, . . . , 9n) 



N-l 

IK- 



CN 



N 

E 



Evaluating the integral in (3.6) as the sum of residues we arrive at the right 
side of dO. □ 



3.2. Schur functions. Here we specialize the formulas of Section 3.1 to the 
Schur functions. 

Proposition 3.4. Rational Schur functions sx{xi, . . . ,xn) (os above we 
identify A G GT n with ^ = X + 5 £ GT n ) are class of determinantal func- 
tions with 



J-i 



g{x; m) = j;"^, a{x) 



1 



1 



/3(x) 



1, 



CAT 



n 



l<i<j<N 



N-l 



fjqx) - f{x) 

Q-l 



Proof. This immediately follows from the definition of Schur functions (2.1 ) 



and evaluation formula (2.2) 



□ 



Propositions 3.2 and |3.3| specialize to the following theorems. 
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Theorem 3.5. For any signature A E GTjv and any k < N we have 

ni=i 11^=1 [Xi-q^ ^) 



det ^ 5a (x.; N, q) U^jx. - g^'^) 

A{xi,...,xk) y [A^-l]^' 



i=l 

where Di^q is the difference operator acting on the function f{xi) by the 
formula 

Theorem 3.6. For any signature A G GTn and any x G C other than or 
q\ i e {0, . . . , N — 2} we have 

where the contour C includes the poles at Xi + N — 1, . . . , Aat and no other 
poles of the integrand. 



Remark. There is an alternative derivation of Theorem 3.6 suggested by 
A. Okounkov. Let x = q^ with k > N. The definition of Schur polynomials 
implies the following symmetry for any fi,X G GT^v 

Using this symmetry, 

5aW ,^V,g - ^^3— , 

where hf^ = 5(^ 0,...) is the complete homogeneous symmetric function. In- 
tegral representation for h^ can be obtained using their generating function 
(see e.g. (Ml Chapter I, Section 2]) 



Extracting kg as 

hp 



oo N ^ 

H{z) = ^he{yi,...,yN)z^ = 11?^ • 

- - ■ ^ ^ zyi 



2^ / 



we arrive at the integral representation equivalent to Theorem 3.6 In 
fact the symmetry (3.8) holds in a greater generality, namely, one can re- 
place Schur functions with Macdonald polynomials, which are their (g, t)~ 
deformation, see [M, Chapter VI]. This means that, perhaps, Theorem |3.6| 
can be extended to the Macdonald polynomials. On the other hand, we do 
not know whether a simple analogue of Theorem 3.5 for Macdonald polyno- 
mials exists. 



Sending g — )• 1 in Theorems 3.5, 3.6 we get 
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Theorem 3.7. For any signature A E GTjv and any k < N we have 

k 



Sx{xi, . . . iV, 1) = 



{N-i)l 



L(iV-l)!(x,-l) 

k 



N-k 



det 



D: 



where Di i is i/ie differential operator x - ^ 



dxi ■ 



Theorem 3.8. For any signature A G GT^v a^^f^ any x € C other than or 
1 we /laue 

(Af-1)! 1 



(3.9) Sx{x;N,l) = -^f- — : (h dz, 

where the contour C includes all the poles of the integrand. 

Let us state and prove several corollaries of Theorem |3.7[ 

For any integers j, i, N, such that < i < j < N , define the polynomial 

Pj/.N{x) as 

(3.10) 



Pj,e,N{x) 



j-l\N\N-j)\ 



{x - 1) 



dx) 



(x-l) 



(7V-1)! 

it is easy to see (e.g. by induction on j — £) that Pj/^n is a polynomial in 
X of degree j — i — 1 and its coefficients are bounded as — )• oo. Also, 
Pj,o,N{x)=xi-^ + Oil/N). 

Proposition 3.9. For any signature A G GT^r and any k < N we have 

(-l)Q) 



(3.11) Sx{xi,...,Xk;N,l) 



A(xi, ...,Xk) 



X det 



'^^DfJSx(xi;N,l)] . 



J jj=l 



Proof. We apply Theorem 3.7 and, noting that 



xiy [f{x)9{x)] = Q (^^)' [/(-)] (-^ 



dx 



[9{x)] 



for any / and we obtain 



(3.12) Sx{xi,...,Xk]N,l) 



1 



A(xi, . . . ,xfc) 



det 



{N-j)\Dl\Sx{xi-N,l){xi-l) 



(iV-1)! 



det 



nZlDU[Sx{xf,N,l)]{ 



{xi - 1)^-'= 



A(xi, ...,xk) 



□ 
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Corollary 3.10. Suppose that the sequence X{N) G GT^v is such that 

lim SxiN){x;N,l) = ^{x) 

A*— s-oo 

uniformly on compact subsets of some region M C C, then for any k 
lim Sx(N){xi, ...,Xk;N,l) = ^>(xi) • • • <^{xk) 

uniformly on compact subsets of . 

Proof. Since S\(^jy[-^{x] N, 1) is a polynomial, it is an analytic function. There- 
fore, the uniform convergence implies that the limit <I>(x) is analytic and all 
derivatives of Sx(N)ix) converge to the derivatives of 



Now suppose that all Xi are distinct. Then we can use Proposition |3.9 

and get as — oo 



det 



SxiN){xi,. . . ,Xk;N,l) 



(Xi - l)''-lSxiN)ix^■, N, l)Pj,0,NiXi) + 0{1/N) 



A(a;i, ...,Xk) 



det 



(x. - l)''-^SxiN)ix^■, N, l)xl^ + 0{1/N) 



A{xi, ...,Xk) 



det 



llSxiN){xi;N,l) 



{xi - If-ix^^ 



+ 0{l/N) 



A(xi, . . . 



\{Sx^N)[x^■N^) (1 + 



i=l 



0{l/N) 



A(xi, ...,Xk) 

where 0{1/N) is uniform over compact subsets of M^. We conclude that 
(3.13) lim Sx(N){xu.--,Xk;N,l) = <^{xi)---<^{xk) 

uniformly on compact subsets of 

\ \]{x, = xj}. 



i<j 



Since the left-hand side of (3.13) is analytic with only possible singularities 

Xi 
□ 



at for all A^, the uniform convergence in (3.13) also holds when some of Xi 
coincide. 



Corollary 3.11. Suppose that the sequence X{N) £ GTjy is such that 

A->oo N ^ ' 

uniformly on compact subsets of some region M C C, in particular, there is 
a well defined branch of logarithm in M for large enough N. Then for any 
k 

In {Sx{N)ixi, ...,Xk;N, 1)) 



lim 

A-s>oo 



N 



^ixi) + --- + '^{xk) 



um 



iformly on compact subsets of M^. 
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Proof. Notice that 

\ dx 



^)^SA(x;iV,l) 



— HSx{x- AT, 1)), . . . , _ ln(5A(x; N, 1)) 



Sx{x;N,l) 

i.e. it is a polynomial in the derivatives of ln(S'A(a;; 1)) of degree j and so 



Sx{x-N,l) 
Thus, when 



exists, then 



converges and so does 



X, \n{Sx{x- N,l)),...,^ HSxix; N, 1)) 



lim 



In (5Aw(x;iV,l)) 



iV 



det 



ni=i 5'A(Ar)(a;i; N, l)A(a;i, ...,Xk) 

Now taking logarithms on both sides of ( |3.11 ), dividing by A^, factoring out 
Y\i=i S\{N){xi; N, 1) and sending — )• oo we get the statement. □ 

Corollary 3.12. Suppose that for some number A 



uniformly on compact subsets of domain D C C as N ^ oo. Then 
(3.14) 



lirn^ Sx(N) (e ^ , . . . , e ^ ; iV, l) exp (^aVn (yi + • • • + y^)) = J] ^ivi) 



i=l 



uniformly on compact subsets of] 



Proof Let 5A(7v)(e^/^; l)e^^y = Gjv(y)- Since G7v(y) are entire func- 
tions, G(y) is analytic on D. Notice that 



X — / (yNln{x)^ = VNf (^ViVln(x)) 



therefore 
d 



xg) S,,,,ix;N,l) = N^/^ 



y=\/N InQ 



r=0 ^ ^ 



y=\/N In X 

e-^^yGN{y) {l + 0{l/^) 



y=\ TV In x 
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since the derivatives of G]y{y) are uniformly bounded on compact subsets 
of D as — )• oo. Further, 

(x - lY = N-^/^y\l + 0{1/Vn)), X = e^/^, 

and Pj-£,Ar (ef/^) = 1 + 0{l/^/N) with 0{l/^/N) uniformly bounded on 



compact sets. Thus, setting Xi = e^*/^ in Proposition 3.9, we get (for 
distinct yi) 



S 



\(N) 



\xi - lf~^GN{y^) (l + 0(1/^^)) 



v(xi, . . . 



1 fc 



det 



GNivi) ■ ■ - GNiyk)- 



k 



{xi - 1^-^ [1 + 0(1/ ViV) 

A(2;i, ...,Xk) 
= GNiyi) ■ ■ ■ GNiyk) (l + Oil/VN)) ■ 

Since the convergence is uniform, it also holds without the assumption that 
yi are distinct. □ 

3.3. Symplectic characters. In this section we specialize the formulas of 
Section |3.1| to the characters xx of the symplectic group. 

For fi £ CT J let 



det 



H-i+l — /X, — 1 



N 

«J=1 



A(xi,...,XAr) 



Al{xi,.. .,Xn) 
Clearly, for A G GT^ we have 

A,^s = Xxixu . . . , x^) ^TTTT^^^p ' 

where xx is a character of the symplectic group Sp{2N). 

Proposition 3.13. Family ^^(xi, . . . ,xn) forms a class of determinantal 
functions with 

„x+l ^—x—l 

m+1 — m— 1 



q\ g{x;m) =x"'+' -x-'^-^ /3(x) 



1 



a{x) 



qx+l ^ _ 2 



j{x+l)/2 _ q-{x+l)/2 



q-1 



(9-1) 



Ar2 



i<i<J<^ ^ (-l)(2)A(g,...,g 



, [Tf]{x) 



fjqx) + f{q~ 
(9-1)2 



Proof. Immediately follows from the definitions and identity 



At{q,...,q" 



i-l)i'2)A{q,...,q-) 



□ 
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Let us now specialize Proposition 3.2 
We have that 

/ N—k\ 

^ I ^ Xx{xi,...,Xk,q,...,q" 
Xa(xi, . . . ,Xk;N,q) = 

^ A,(g, . . . , q^)A{xi, ...,Xk,q,..., g^-^) A^^jx^, . . . , x^, q, . . . , g^"^) 
A,{x^,...,Xk,q..., q^->')A{q, . . . , g^) Af,{q, ...,q^) 

Theorem 3.14. For any signature X € GT^ and any k < N we have 



A, {q,...,q^){q-lf-H -lP 
As{xi, . . . ,Xk,q . . . ,q'' 



(3.15) XA(xi,...,x,;jV,g) = --W'--''^ y 



i=l 

where j is the difference operator 

f{qx)+f{q-'x)-2f{x 



fix) 



acting on variable Xi. 



Remark. Note that in Proposition 3.13 the difference operator differed 
by the shift 2/(g — 1)^. However, this does not matter, as in the end we use 
the operator \\i^;{Dli - Dlj). 

Proposition |3.3| after algebraic manipulations to compute the coefficient 
in front of the integral yields. 

Theorem 3.15. For any signature X G GT^ and any q ^ 1 we have 
(3.16) Xa)- (-1>"-''°W('-1)="-'PA'1,1 



{xq;q)N-iix ^q; q)N-i{x - x ^)[N]q 

1 r (x^+i - x-^-1) 

27ri / ^ 

i=l 



-dz 



with contour C enclosing the singularities of the integrand at z = Xi+ N 
1, . . . , Aat. 



Theorem 3.15 looks very similar to the integral representation for Schur 



polynomials, this is summarized in the following statement. 
Proposition 3.16. Let X G GT^, we have 

Xxix;N,q) = ^^^5,(xg^-i;2iV,g), 
X + 1 

where v G GT2Ar is a signature of size 2N given by Vi = Xi + 1 for i 
1, . . . ,N and Vi = —X2N-i+i for i = + 1, . . . , 2N . 

Proof. First notice that for any /i G GT^, we have 



(x^ - 



ic 



liM'' + - g~^'~^ - g^'+^) Ic Yiiia' + - q^^'~^ 
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where C encloses the singularities of the integrand atz = Ai + A^ — l,...,AAr 
and C encloses all the singularities. Indeed, to prove this just write both 
integrals as the sums or residues. Further, 



+ q ^ 



Therefore, the integrand in (3. 16) transforms into 
(3.17) 



X q 



Nz 



_ qX,+N-i+i^(_^qz _ ^-(A,+Af-i)-i) 



z' z' -N 
a X q 



where z' = z + N. T he contour integral of (3.17) is readily identified with 
that of Theorem 3.6 for S,^{xq^^^;2N,q). It remains only to match the 
prefactors. □ 

Next, sending q ^ I we arrive at the following 3 statements. 
Define 



(3.18) Ai(xi,...,Xfc,l^-'=) = lim 



Asjxi, ...,Xk,q,...,q^ ^) 
(g-l)( 2 ) 



As(xi, . . . , X,) n JT"^ n - - k)i 



Theorem 3.17. For any signature X E GT^ and any k < N we have 



l<i<j<N-k 
N 

AKl^) 



(3.19) X.(xi,...,x,;jV,i) = ^^^^^^^^^^-^^.^^_^^ x 



(-l)^2;det 



'' dx, 



^ k 

llXx{xi;N,l) 

i,j=i j=i 



(xj-Xj ^)i2-Xi-X i 1) 
2(2A^- 1)! 



l\Af-l 



Remark. The statement of Theorem 3.17 was also proved by De Gier 
and Ponsaing, see jGPj . 

Theorem 3.18. For any signature X G GTj!^ we have 



Xa(x;A^,1) 



2(2iV- 1)! 



{x - x-^){x + x-^ -2)^-1 



1 



[x^ — X ^) 



-dz, 



27ri Jc UZii^ -{Xi + N-i + l)){z + Xi + N-i + l)' 
where the contour includes only the poles at Xi + N — i + 1 for i = 1, . . . , N . 
Proposition 3.19. For any signature X G GT^ we have 



(3.20) 



:^A(x;iV,l) 



-Suix;2N,l), 



x + l 

where v E GT2Ar is a signature of size 2N given hy Vi = Xi + 1 for i 
l,...,N andvi = -X2N-i+i for i = N + 1, . . . , 2N. 
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Remark. We believe that the statement of Proposition 3.19 should be 
known, but we are unable to locate it in the literature. 

Analogously to the treatment of the multivariate Schur case we can also 



derive the same statements as in Proposition 3.9 and Corollaries 3.10 3.11 



3.12 for the multivariate normalized symplectic characters. 



3.4. Jacobi polynomials. Here we specialize the formulas of Section 3.1 



to the multivariate Jacobi polynomials. We do not present the formula for 



the version of (|2.5|), although it can be obtained in a similar way. 



Recall that for A G GT+ 



J\{zi, . . .,Zk;N,a,b) 



^xjzi, ...,Zk,l 

3A(l^;a,6) 



We produce the formulas in terms of polynomials /i G GT^ and, thus, 
introduce their normalizations as 



P^{xi,. . . ,Xk;N,a,b) 



q}„ (xi,...,Xk,l 



N-k. 



These normalized polynomials are related to the normalized Jacobi via 



J\{zi, ...,Zk; N,a,b) = Pf,{ 



z\ + Zi 



-1 



where as usual Aj + 



/ij for i = 1, . . . , A^. 



;Af,a,6), 



Proposition 3.20. The polynomials P^{xi, . . . ,xn), /x G GT^ are a class 
of determinantal functions with 

T{x + a + l) 



1, g{x;m) = pm{x; a,b), a{x) = x{x+a+b+l), I3{x) 



r(x + l)r(a) 



N 



CN 



n r(r)r(a) -|-|- 
T(r + a) J-J- 

r=l ^ ^ l<i<j<N 

92 



{j -i){2N -i- j + a + b + l) 
d 



r= (x2-l)^ + ((a + 6 + 2)x + a-6))„ 
ax^ ox 

Proof. We have (see e.g. [002^ Section 2C] and references therein) 
(3.21) P^(l"; a, 6) = n ^^^^l^i)^^ >< " + l^j + a + b + 1) 



n 



r(r) 



n 



J- r(r + a) 

r=l ^ ' 0<i<j<n 

and also (see e.g. |Edj . |KoSwj ) 

m{m+2a)pmix;a,b) 



1 

{j -i){i + j + a + b+l)'' 



q2 q ' 

Now the statement follows from the definition of polynomials Pn 



Pni{x;a,b), 
□ 



d_ 
dx 



Specializing Proposition 3.2, using the fact that for x = — we have 



1 — Qz 



and Pfi{x) = J\{z), we obtain the following. 
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Theorem 3.21. For any A G GT^ and any k < N we have 
(3.22) Jx{zi,...,Zk;N,a,b) 

_ r(m + a)r(2m- 1 + + 6) 1 

det[Pl;a,=i ^ ^ ^^ (z, + zri-2)^"^r(iV + a + 6) 

2(^) A(z, + . . . , + ^,-) M ' r(A^ + a)r(2iV - 1 + a + 5) 

where T)i^a,b is the differential operator 

2 9^ {{a + b + 2){zi + zr^) + 2a-2b-2z-^) 8 



Next, we specialize Proposition 3.3 to the case of multivariate Jacobi 



polynomials. Note that thanks to the symmetry under ^ + (a + 6 + l)/2-f-)- 
— (C + (a + ft + 1) /2) of the integrand we can extend the contour C to include 
all the poles. 

Theorem 3.22. For any A G GT^ we have 
(3.23) 

T ^ AT ^^ 5A(x,l^-i;a,6) r(2iV + a + 6-l 1 
Jx{z;N,a,b) ^ - 



^x{l^; a, b) r(n + a + b)T{a + 1) /^+^ _ 



2 

2^ 



-C,C + a + 6 + l;a + l;-^^^^ )(C + (a + 6 + l)/2) 



2vri/c ni(C-/^i)(C + Mi + a + fe+l) 

where the contour includes the poles of the integrand at C, = — (a + 6+l)/2ib 
{Hi + {a + b+ l)/2) and Hi = Xi + N - i for i = 1, . . . , N. 

4. General asymptotic analysis 

Here we derive the asymptotics for the single-variable normalized Schur 
functions Sx{x] N,l). In what follows O and o mean uniform estimates, not 
depending on any parameters and const stays for a positive constant which 
might be different from line to line. 

4.1. Steepest descent. Suppose that we are given a sequence of sig- 
natures A(iV) G GTjv (or, even, more generally, X{Nk) G GTat^. with 
Ni < N2 < N3 < . . .). We are going to study the asymptotic behavior of 
S\(N)ix] A^, 1) as — )• cxD under the assumption that there exists a function 
f{t) for which as — ^ 00 the vector (Ai(A^)/A^, . . . , Xn{N)/N) converges to 
{f{l/N), . . . , f{N/N)) in a certain sense which is explained below. 

Let Ri , Roo denote the corresponding norms of the difference of the vectors 
(A,(iV)/iV) and fij/N)): 

N 

Ri{XJ) = Y, 



N 



-Roo(A,/)= sup 

i=l...,Af 



urn) 



N 
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In order to state our results we introduce for any y £ C the equation 

^'■'^ Jo = 



Note that a solution to (4.1) can be interpreted as an inverse Hilbert trans- 



form. We also introduce the function T{vu; f) 

(4.2) F{w; f) = C \n{w - f{t) - 1 + t)dt. 

Jo 



Observe that (4.1) can be rewritten as T'{w; f) = y. 



Proposition 4.1. Fory £ M\{0}, suppose that f{t) is piecewise- continuous, 
Roo{X{N), f) is bounded, Ri{X{N),f)/N tends to zero as N ^ oo, and 
Wo = wo{y) is the (unique) real root of ( |4.1[ ). Let further y G M \ {0} be 
such that Wo is outside the interval [ ^^j^^ + 1] for all N large enough. 

Then 

ln5'A(Ar)(e2'; A^, 1) ywo-F{wo) 

lim ^ — = ^ — . 

Af^oo N e(ef - 1) 

Remark 1. Note that piecewise-continuity of f{t) is a reasonable as- 
sumption since / is monotonous. 

Remark 2. A somehow similar statement was proven by Guionnet and 
Maida, see [UMl Theorem 1.2]. 



When f(t) is smooth, Proposition 4.1 can be further refined. For w £ C 
denote 

Proposition 4.2. Let y e R \ {0} be such that wo = wo{y) (which is the 
(unique) real root of (4.1 )j is outside the interval [^^^r/^, ^4?^ + 1] for all 



large enough N . Suppose that for a twice- diff'erentiable function f{t) 
(4.3) hm exp {Q{w; \{N), /)) = g{w) 

uniformly on an open A4 set in C, containing wo. Assume also that g{wo) ^ 
and T" [wo] f) 7^ 0. Then as N ^ 00 



^ (y.r,, ^0 - /(O) - 1 , exp(jV(^^;o-.F(^;o;/))) ^^ , 
'^''^ = V--F"(^o;/)(-o-/(l))^^"°^ e-iey - 1)-- ['^' 

The remainder o{l) is uniform overy belonging to compact subsets o/R\{0} 
and such that wo = WQ{y) £ Ai. 

Remark. If the complete asymptotic expansion of exp{Q{w; X{N), /))) 
as — >■ 00 is known, then, with some further work, we can obtain the 
expansion of Sx(n){^''^'i up to arbitrary precision. In such expansion. 



0(1) in Proposition 4.2 is replaced by a power series in N^^^'^ with coefficients 
being the analytic functions of y. The general procedure is as follows: we 
use the expansion of ex.p{Q{w, X{N), f))) (instead of only the first term) 
everywhere in the below proof and further obtain the asymptotic expansion 
for each term independently through the steepest descent method. This 
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level of details is enough for our applications and we will not discuss it 
any further; all the technical details can be found in any of the classical 
treatments of the steepest descent method, see e.g. [Coj . |Erl] . 

Proposition 4.3. Suppose that f{t) is piecewise-dijjerentiable, 
RociX{N),f) = 0(1) (i.e. it is bounded), and Ri{\{N), f)/VN goes 
to as iV —7- oo . Then for any fixed /i G M 

5A(7V)(e"/^; iV, 1) = exp (^VNEif)h + ^Sif)h'' + o(l)) 

as N ^ oo, where 

Jo Jo Jo 

Moreover, the remainder o{l) is uniform over h belonging to compact subsets 
o/M\0. 

We prove the above three propositions simultaneously. First observe that 
we can assume without loss of generality that Xi^{N) > 0. Indeed, when 
we add some positive integer M to all coordinates of X{N), the func tion 
5>.(- j v)(e^ ; A^, 1) is multiplied by e*^^, but the right sides in Propositions 
|4.2[ |4.3| also change accordingly. 

We start investigating the asymptotic behavior of the integral in the right 



4.1 



side of the integral representation of Theorem 3.8 



(4.4) SA(e^;A^,l) = , A?' , 1 -^j dz 

(e^ - 1)^-' 2vri 7^7 nf=i(^ - (A,(iV) + - i)) 



Changing the the variables z = Nvu transforms (4.4) into 



(4 5) A^[^}}Ln^~n± I exp(iV,^.) 

^ ' (e^-ir-^ 2vriJcnf^i(^-^-^%+^-^) 

From now on we study the integral 
(4.6) / eMNyw) 



where the contour C encloses all the poles of the integrand. 
Write the integrand as 

/ 

exp Nyw - ^ ln(w; - (A^- + N - j) /N) 

Unless otherwise stated, we choose the principal branch of logarithm with a 
cut along the negative real axis. 
Observe that 

Q{w; A, /) = In U; - ^ - ^ Hw - f{j/N) - 1 + j/N) 

i=i ^ ^ i=i 
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Under the conditions of Proposition 4.1 and basing on the approximation 
ln(l + a) w Q we have 



\Q{w;X,f)\=o{N) sup 



w — a 



where A is the smallest interval in M containing the four points 
inf suptg[o,i] (fit) - 1 + infte[o,i] (/(*) - 1 + 0). minj 



Under the conditions of Proposition |4.3| we have 

1 



\Q{w;X,f)\=o{VN) sup 



w — a 



Using a usual second-order approximation of the integral (trapezoid for- 
mula) we can write 



N-l 



N 



^ Hw - f{j/N) + (1 - j/N)) = N\Y, 



ln{w- f{j/N)-l + j/N) 



j=0 



N 



N 



ln{w - f{t) - 1 + t)dt + 



Hw- f{0)-l)-ln{w- f{l)) 



and depending on the smoothness of / we have the following estimates on 
the remainder T: 

(1) If / is piecewise-continuous, then 

T{w, f, N) = o{N) sup I ln{w - f{t) -l + t)\, 



(2) If / is piecewise-differentiable, then 

w-fit)-l + t 



T{wJ,N) = 0{l) sup 

\t-s\<l/N 



In 



w — f{s) — 1 + s 



+0(1) sup 



/'W + 1 



w-f{t)-l + t 



(3) If / is twice differentiable (on the whole interval), then 



T{wJ,N) = 0{l/N)snp 

t 



fit) + 1 



d 



dt \w- f{t) -t + t 



The integral (4.6) transforms into the form 
(4.7) 



exp yN{yw — F{w; f))j r]\i{w)dw. 
Here 

, , /ln(ii;-/(0)-l)-ln(u;-/(l)) , ^, , 

r^iw) = exp ^ '-^ ^ ^ + Q{w, X, f) + T{w, /, N) 



has subexponential growth as — )• oo. 

Note that Re {F{w, /)) is a continuous function in w, while Im {T{w; /)) 
has discontinuities along the real axis, both these functions are harmonic 
outside the real axis. 

The asymptotic analysis of the integrals of the kind (4.7) is usually per- 
formed using the so-called steepest descent method (see e.g. [Coj . |Erl| .). We 
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will deform the contour to pass through the critical point of yw — F{w] /). 
This point satisfies the equation 



(4.8) Q = {yw-F{wJ))' = y 



dt 



w-f{t)-l + t 



In general, equation ( |4.8| (which is the same as (4.1 )) may have several roots 
and one has to be careful to choose the needed one. 

In the present section y is a real number which simplifies the matter. If 



y > then (4.8) has a unique real root WQ{y), moreover wo{y) > 0. Indeed 



the integral in (4.8) is a monotonous function oiw& M+ changing from +co 
down to zero. In the same way if y < then there is a unique real root 
wo{y) and it satisfies wo{y) < 0. Moreover, WQ^y) — )■ oo as y — )• 0. In what 
follows, without loss of generality, we assume that y > 0. 

Next, we want to prove that one can deform the contour C into C which 
passes through wo{y) in such a way that Re{yw — F{w; /)) has maximum 
at WQ{y). When y is real, the vertical line is the right contour. Indeed, (for 
> 2) the integrand decays like as s — )• oo on the vertical line and 
exponentially decays as Rew — t- —oo, therefore the integral over this vertical 
line is well defined and the conditions on wq guarantee, that the value of 
the integral does not change when we deform the contour. Moreover, if 
w = WQ + is, s G M, s ^ then immediately from the definitions follows that 

Re{yw - T{w; /)) < Re{ywo - T{wo; /)). 

Now the integrand is exponentially small (compared to the value at wo) 
everywhere on the contour C' outside arbitrary neighborhood of wo- Inside 
the neighborhood we can do the Taylor expansion for yw — J-{w; /). 

Denoting u = —i^y—J^"{wo; /) and w = wo + s/{u\/lV), the integral turns 
into 



(4.9) exp\^N{ywo-J^iwo;f)) 

X / expi-NT"{wo;f){w-wof/2 + N5{w-wof]rNi'w)dw 

JwQ—ie \ / 

exp ( N{ywo - J^{wo; /)) 



uVN 

^e\u\ 



exp (-sV2 + s^d/y/N) rN{wo{y) + s/{uVN))ds 

-^/Ne\u\ ^ ' 



2vr^^rAr(u;o) exp [N{ywQ - F{wq; /))) , 
uv N 



where 

\6\< sup F"'{w;f) 

and 

l^l = l'^IH-'- 

When we approximate the integral over vertical line by the integral over 



the e- neighborhood (reduction of (4.7) to the first line in (4.9)) the relative 
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error can be bounded as 
(4.10) 

constxexp{N Re{J'{wo+ie; f)—J'{wo; /)) w constxexp{—Ne^\J'"{wo; /)|/2) 



Next, we estimate the relative error in the approximation in (4.9) (i.e. the 
sign in ( |4.9| )). Suppose that e < \u/6\/2 and divide the integration seg- 
ment into a smaller subsegment \s\ < N~^/^^ \J~\/WJ[^\ and its complement. 

When we omit the term in the exponent we get the relative error at most 
const X iV~^/^^ when integrating over the smaller subsegment (which comes 

from the factor exp (^sH/^fN^ itself) and const x exp (^—N^'^/^^\5\^'^/'^/^ 

when integrating over its complement (which comes from the estimate of 
the integral on this complement). 

When we replace the integral over [— -v/iVe|'u|, +\/iVe|?^|] by the integral 
over (— oo,-|-oo) in (4.9) we get the error 

const exp ( — Ne^ \u^\/2), 

Finally, there is an error of 

const sup \rN{w) — r]\j{wQ)\ 

coming from the factor r]\f{w). Summing up, the total relative error in the 
approximation in (4.9) is at most constant times 



(4.11) Ar-3A0 + exp f-iV-2/i5|5|-2/3/4\ + e^p(-Ne^\u^\/2 



+ sup \rN{w) - rN{wo)\. 
w£[wo—ie,wo+ie] 



Combining (4.5) and (4.9) we get 

gA(eM^-^) 
5a(1^) 
1 (A^-1)! 



27r {ey - 1 



N 



l-N 



Using Stirling's formula we arrive at 
(4.12) 



rN[wo) 



«a(1^) 



rNiwo) exp{N{ywo-J^iwo; /))). 



exp{N{ywo - J^iwo; /))). 



With the relative error in ( |4.12 ) being the sum of (4.10), (4.11) and 0{1/N) 
coming from Stirling's approximation, and £ satisfying e < \u/5\/2. 

Now we are ready to prove the three statements describing the asymptotic 
behavior of normalized Schur polynomials. 



Proof of Proposition 4-1 Use (4.12) and note that after taking logarithms 
and dividing by N the relative error in (4.12) vanishes. □ 



Proof of Proposition ^.2. Again this follows from (4.12). It remains to check 
that the error term in (4.12) is negligible. Indeed, all the derivatives of 
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as well as \5\, \5\ are bounded in this limit regime. Thus, choosing 
e = N^^^^^ we conclude that all the error terms vanish. □ 



Proof of Proposition 4-3. The equation (4.8) for wq reads 

0. 



Clearly, as — t- 
dt 



h/VN 
oo we have wq 



wo-f{t)-l + t 

N /h ^ CO. Thus, we can write 



Wo 



Denote 




1 + t 



Wo 



fit) + l-t 



Wo 



+ 0{l/{wof) dt. 



A= [\f{t) + l-t)dt, B= f\f{t) + l-tfdt 
Jo Jo 



and rewrite (4.8) as 



Wo 



wo- 



h 



0(1), 



If follows that as — 7- oo we have 



wo = VN/{2h) + 



^+4^ + 0(l/V]V) 



/l2 



h 



h 



+ A + 0{1/VN). 



Next, let us show that the error in (4.12) is negligible. For that, choose e 
in ([49j) to be N^/^°. Note that \ J^"{wo;f)\ is of order A^-^, and \T"'{w;f)\ 
(and, thus, also \6\) is of order N~^/'^ on the integration contour and \u\ is of 
order iV"^/^. The inequality e < \u/5\/2 is satisfied. The term coming from 
(4.10) is bounded by exp(— const x A^^/^) and is negligible. As for ( |4.11[ ) 
the first term in it is negligible, the second one is bounded by exp(— const x 
N"^/^^) and negligible, the third one is bounded by exp(— const x N^l^) 
which is again negligible. Turning to the fourth term, the definition of 
r]\f{w) implies that both r^iw) and r^iwo) can be approximated as 1 + o(l) 
as A^ — )• oo and we are done. 

Note that 



1 



^/-r'iwo; f) 



l + o(l) 



as A" — ^ oo. Now (4.12) yields that 
(4.13) 

SA(e^/v^, i^-i; 

As A^ — )• oo we have 



exp 



N{-l-\n{e''l^'' -l)+hwo/VN-F{wo; /)) (l+o(l)) 



^{wo; f) 

= ln(u;o) + 



f \n{wQ - f{t) - 1 + t)dt 
Jo 

fit) (fit) + 1 - t)2 



t- 1 



Wo 



2wl 



+ 0{l/w, 



dt 
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In(ifo) 



A 
Wo 



B 

2^ 



+ 



Wn 



and using (4.8) 



hwo 



1 + 



A 

Wo 



+ 



B 



+ 



Wn 



Thus, 



1 - ln(e''/^ - 1) + h/^/NwQ - T{wo; f) 



Wq 



w, 







Wo 



2wl 



In [woh/\fN 



In 1 + 



A B ( A B 

2w^ \wq w^ 



Wq 



2VN 
h 



+ — + 0{N 3/2) + + . 



6A^ 

/l2 



A B/2 + A^ 

-+ 2 



1 / h 

m^2 V 2~ 



two ' 2wl 



-0{N 



-3/2N 



Wo 



W, 







/l2 



2~'n 



h 



2,/N 24N 



+ 0(iV-3/2) 



To finish the proof observe that 



A = E{f) + 1/2, B= [' f{t)dt + 2 C f{t){l - t)dt + 1/3, 

Jo Jo 

thus, (4.13) transforms into 

eME{f)hVN + S{f)/2)il + o{l)). 



□ 



4.2. Values at complex points. The propositions of the previous section 
deal with S\{e'^; N,l) when y is real. In this section we show that under 
mild assumptions the results extend to complex ys. 

In the notations of the previous section, suppose that we are given a 
weakly-decreasing non-negative function f{t), the complex function J-'{w; /) 
is defined through ( |4.2[ ) , y is an arbitrary complex number and wq is a critical 
point of yw — T{w; /), i.e. a solution of equation ( |4.8| ). 

We call a simple piecewise-smooth contour 7(5) in C a steepest descent 
contour for the above data if the following conditions are satisfied. 

(1) j{0)=wo, 

(2) The vector {F"{wq; f))~^^'^ is tangent to 7 at point 0, 

(3) Re{y'y{s) — T{'y{s); /)) has a global maximum at s = 0, 

(4) The following integral is finite 

j exp (i?e(y7(s) - -^(7(5);/)) h'it)\dt < 00. 

Remark. Often the steepest descent contour can be found as a level line 
Im{yw - T{w; /)) = Im{wo - T{wo; /)). 
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Example 1. Suppose that f{t) = 0. Then 

J^{w; f)= [ ln{w - 1 + t)dt = w ln{w) - {w - 1) ln{w - 1) - 1 

and 

(4.14) T'{w; 0) = ln(u;) - ln(u; - 1) = - ln(l - l/w). 

And for any y such that 7^ 1, the critical point is 

Wo = wo{y) = 

Suppose that is not a negative real number, this implies that wq does 
not belong to the segment [0, 1]. Figure [7] sketches the level lines Re{yw — 
J^{w;0)) = Re{ywQ — F{wq;Q)). In order to understand that the picture 
looks like that, observe that there are 4 level lines going out of wq. Level lines 
can not cross, because yw—F{w, 0) has no other critical points except for wq. 
When \w\ S> 1, we have Re{yw — F{w; 0) « Re{yw) — In {wl, therefore level 
lines intersect a circle of big radius Rm2 points and the picture should have 
two infinite branches which are close to the rays of the line Re{yw) = const 
and one loop. Since the only points where ex.p{Re{y'w — J^{w;0)) is not 
analytic form the segment [0, 1], the loop should enclose some points of this 
segment (real part of analytic complex functions without critical points can 
not have closed level lines). 




Figure 7. Sketch of the level lines Re{yw - T{w;0)) = 
Re{ywo - J^{wq; 0)) for y = 1 - 1. 

Now the plane is divided into three regions A, B and C as shown in 
Figure [?} Re{yw — T{w;0) > Re{ywo — T{'Wo;0) in A and C, while 
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Re{yw — F{w,Q) < Re{ywo — J^(tt;o;0), this can be seen by analyz- 
ing Re(yw — J-{w,0) for very large \w\. There are two smooth curves 
Im{yw — T{w; 0) = Im{ywo — F{wq] 0) passing through wq. One of them 
has a tangent vector parallel to ^^T"{wq; 0) and another one parallel to 
i^J^^^Jwq]^ . Take the former one, then it should lie inside the region B. 
In the neighborhood of wq this is our steepest descent contour. The only 
property which still might not hold is number 4. But in this case, we can 
modify the contour outside the neighborhood of wq, so that Re{yw — T{w; 0) 
rapidly decays along it. Again, this is always possible because for \w\ 3> 1, 
we have Re{yw — F{w, 0)) w Re{yw) — 1 — In |tt;|. 

Example 2. More generally let f{t) = a{l — t), then 

F{w]a{l-t))= ! \n{w + {a + l){t - l))dt 
Jo 

w ln{w) — {w — {a + 1)) ln(u; — (a + 1)) ^ 
~ a + 1 

and 

ln(l - {a + l)/w) 



T'{w;a(l - t)) 



a + 1 

For any y such that ^ 1, the critical point id 

wo = woiy) = ia + l)/il-e-y^''-^^^). 

Note that if we set w = u(a + 1), then 

T{w; a{l - t)) = uln{u) - {u - 1) ln(n - 1) + ln(a + !)-!, 

which is a constant plus J-{u;0) from Example 1. Therefore, the linear 
transformation of the steepest descent contour of Example 1 gives a steepest 
descent contour for Example 2. 

Proposition 4.4. Suppose that f{t), y and wq are such that there exists 
a steepest descent contour 7 and, moreover, the contour of integration in 
(4.4) can be deformed to 7 without changing the value of the integral. Then 
Propositions 4-1 and 4-2 hold for this f{t), y and vuq. 



Proof. The proof of Propositions 4.1 and |4.2| remains almost the same. The 
only changes are in formula (4.9) and subsequent estimates of errors. Note 
that condition 4 in the definition of steepest descent contour guarantees that 
the integral over 7 outside arbitrary neighborhood of wq is still negligible as 
00. 

Observe that the integration in (4.9) now goes not over the segment [wq — 
ie, wo+ie] but over the neighborhood of wq on the curve 70. This means that 
in the relative error calculation a new term appears, which is a difference of 
the integral 

'^y^ds 



over the interval [— v Ae|M|, — v Ae|ii|] of real line and over the part of 



rescaled curve inside circle of radius VNelul around the origin. 

The difference of the two integrals equals to the integral of exp(— over 
the lines connecting their endpoints. But since 1/u = —{F"{wQ]f))~^^'^ 
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is tangent to 7 at 0, it follows that for small e the error is the inte- 
gral of exp(— over segment joining \/iVe|ti| and ^/~Ne\u\ + Qi plus 
the integral of exp(— joining —y/Ne\u\ and — \/iVe|ii| + Q2 with 
IQil < (\/iVe|u|)/100 and similarly for Q2. Clearly, these integrals expo- 
nentially decay as — )• 00 and we are done. □ 



It turns out that in the context of Proposition 4.3 the required contour 
always exists. 

Proposition 4.5. Proposition\4.3\ is valid for any h £ C. 



Proof. Recall that in the context of Proposition 4.3 y = h/\/N and goes 
to as A'' — )• 00, while wq k, 1/y goes to infinity. In what follows without 
loss of generality we assume that h is not an element of M<o. (In order to 
work with h G M<o we should choose other branches of logarithms in all 
arguments.) 

Let us construct the right contour passing through the point wq. Choose 
positive number r such that r > \ f{t) + 1 — t\ for all < t < 1. Set ^ to be 
the minimal strip (which is a region between two parallel lines) in complex 
plane parallel to the vector \/h and containing the disk of radius r around 
the origin . 

Since wq is a saddle point of yw — F{w\ /), in the neighborhood of wq the 
are two smooth curves Im{yw — T{w; /)) = Im{ywo — F{wq] /)) intersecting 
at wq. Along one of them Re{yw — F{w; /)) has maximum at wq, along 
another one it has minimum; we need the former one. Define the contour 
7 to be the smooth curve Im{yw — F{w]f)) = Im{ywoT{wo; f)) until it 
leaves and the curve (straight line) Re{yw) = const outside 

Let us prove that Re{yw — F{w; /)) has no local extremum on 7 except 
for wq, which would imply that wq is its global maximum on 7. First note, 
that outside ^ we have 



Re{yw - F{w; /)) = Re{yw) - [ ln\w - f{t) + 1 - t\dt, 

Jo 

with the first term here being a constant, while the second being monotonous 
along the contour. Therefore, outside ^ we can not have local extremum. 
Next, straightforward computation shows that if is large enough, then one 
can always choose two independent of N constants 1/2 > Gi > and G2 > 
such that Re{yw — F{w]f)) > Re{ywQ — F{wo; f)) for w in satisfying 
= Gi|t(;o| or {wl = G2\wo\. It follows, that is Re{yw — F{w; f)) had 
a local extremum, then such extremum would exist at some point wi £ ^ 
satisfying Gi|wo| < l^^il < G2|t«o|- But since Im{yw — F{w, /)) is constant 
on the contour inside we conclude that wi is also a critical point of 
yw — F{w; /). However, there are no critical points other than wq in this 
region. 



Now we use the contour 7 and repeat the argument of Proposition 4.3 



using it. Note that the deformation of the original contour of (4.4) into 



7 does not change the value of the integral. The only part of proof of 



Proposition 4.3 which we should modify is the estimate for the relative error 



in (4.12). Here we closely follow the argument of Proposition 4.4 The 



only change is that the bound on Qi and Q2 is now based on the following 
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observation: The straight hne defined by ?liyw = ?R-ywQ (which is the main 
part of the contour 7) is parallel to the vector i/y. On the other hand, 

VF"{wo) = i/y(l + Oil/VN) « i/y. 

□ 



Remark. In the proof of Proposition 4.5 we have shown, in particular, 
that the steepest descent contour exists and, thus, asymptotic theorem is 
valid for all complex y, which are close enough to 1. This is somehow similar 
to the results of Guionnet and Mai'da, cf. \GM.\ Theorem 1.4]. 

5. Statistical mechanics applications 

5.1. GUE in random tilings models. Consider a tiling of a domain 
drawn on the regular triangular lattice of the kind shown at Figure [T] with 
rhombi of 3 types which are usually called lozenges. The configuration of 
the domain is encoded by the number which is its width and integers 
f^i > > ' ' ' > fJ-N which are the positions of horizontal lozenges sticking 
out of the right boundary. If we write /ij = Xi + N — i, then A is a signa- 
ture of size A^, see left panel of Figure [TJ Due to combinatorial constraints 
the tilings of such domain are in correspondence with tilings of a certain 
polygonal domain, as shown on the right panel of Figure [TJ 

Let Qx denote the domain encoded by A S GTjv and define T;^ to be a 
uniformly random lozenge tiling of 0,x. We are interested in the asymptotic 
properties of T\ as — )• 00 and A changes in a certain regular way. 

Given T\ let vi > 1^2 > ' ' ' > i^k be horizontal lozenges at kth vertical 
line from the left. (Horizontal lozenges are shown in blue in the left panel 
of Figure [ij) We again set Vi = Ki + k — i and denote the resulting random 
signature n of size k via T^. 

Recall that the GUE random matrix ensemble is a probability mea- 
sure on the set of k x k Hermitian matrix with density proportional to 
exp(— Trace(X^)/2). Let GUEfc denote the distribution of k (ordered) eigen- 
values of such random matrix. 

In this section we prove the following theorem. 

Theorem 5.1. Let X{N) G GTjy, A^ = 1, 2, . . . be a sequence of signatures. 
Suppose that there exist a non-constant piecewise-differentiahle weakly de- 
creasing function f{t) such that 

" taiN) =„(ViV) 



E 



A^ 



1=1 

as N ^ 00 and also supj ^ |Aj(A^)/A^| < 00. Then for every k as N ^ 00 
we have 



^iN) - NE{f) 



GUE 



in the sense of weak convergence, where 

E{f) = C fit)dt, Sif) = C f{tfdt - E{ff + f{t){l - 2t)dt. 
Jo Jo Jo 

Remark. For any non-constant weakly decreasing f(t) we have S{f) > 0. 
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Corollary 5.2. Under the same assumptions as in Theorem \5.1\ the 
(rescaled) joint distribution of k{k + l)/2 horizontal lozenges on the left 
k lines weakly converges to the joint distribution of the eigenvalues of the k 
top-left corners of a k x k matrix from GUE ensemble. 

Proof. Indeed, conditionally on the distribution of the remaining k{k — 
l)/2 lozenges is uniform subject to interlacing conditions and the same prop- 
erty holds for the eigenvalues of the corners of GUE random matrix, see [Bar] 
for more details. □ 



Let us start the proof of Theorem 5.1 



Proposition 5.3. The distribution of is given by: 

Prob{T, = ,} = , 

where sxj^ is the skew Schur polynomial. 

Proof. Let k G GTm and fj, G GTj\,/_i. We say that k and yu interlace and 
write /X ^ At, if 

Ki > fl2> K2 > ■ ■ ■ > fJ-M^l > KM. 

We also agree that GTq consists of a single point point, empty signature 
and ^ K for all k E GTi. 

For K G GTk and /x E GT^ with > L let Dim(/i, k) denote the number 
of sequences -< C^~^^ -< ••• -< such that C G GTj, = k and 
= fi. Note that through the identification of each with configuration 
of horizontal lozenges on a vertical line, each such sequence corresponds to a 
lozenge tiling of a certain domain encoded by n and /i, so that, in particular 
the tiling on the left panel of Figure [T] corresponds to the sequence 

0^ (1) -< (3,-1) ^ (2,0,-1) -< (2,2,-1,-1) ^ (3,2,2,-1,-1). 

It follows that 

ProbjT, - r?) - j3im(0. A) " 

On the other hand the combinatorial formula for (skew) Schur polynomials 
(see e.g. [Ml Chapter I, Section 5]) yields that for k £ GTk and /i G GT^ 
with K > L we have 

Dim(/i, k) = ^^(1"^"^)' = ^m(I^)- 

□ 

Introduce multivariate normalized Bessel function Bk{x;y), x = 
(xi, . . . , Xfc), y = {yi,...,yk) through 



detij=i,,..,fcf exp(xiyj) 

^"^''"^ = n.<,fe-^)n.<,te-«) n*^- - 

The functions Bk{x\y) appear naturally as a result of computation of 
Harish-Chandra-Itzykson-Zuber matrix integral (1.11). Their relation to 
Schur polynomials is explained in the following statement. 
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Proposition 5.4. For A = Ai > A2 > • • • > G GTk we have 

'-TT]^ Jl ^ = Bk{xi, ...,Xk;Xi + k-l,X2 + k-2...,Xk) 

S\[L ) .^^ Xi Xj 

Proof. Immediately follows from the definition of Schur polynomials and the 
evaluation of s\{l^) given in (2.3). □ 

We study T| through its moment generating functions KBk{x; + 6^), 
where x = (xi, . . . , x^), 6k = (fc — 1, /c — 2, . . . , 0) as above and E stays for the 
expectation. Note that for k = 1 the function Ei?fc(x;T^ + 5^) is nothing 
else but usual one-dimensional moment generating function Eexp(xT^). 

Proposition 5.5. We have 



sa(I^) Xi 

^<i<j<k 



Proof. Let Z = (zi, . . . , Zm) and Y = (yi, 
(see e.g. [Ml Chapter I, Section 5]) 

s^iZ).s^/,{Y) = s,iZ,Y) 



yn) and let /x G GT m+n, then 



Therefore, Propositions |5.3 and 5.4 yield 

3C j X '■j 



(ES,(x;TU4))n5- 



E 



5-.) .,(l^).,/r,(l^-') 



sa(I^) 



E 



s,(e^\ . . . , e-'=)s;,/,(l^^-'=) SA(e-^ . . . , 1^-^=) 



^a(I^) 



«a(1^) 



□ 



The counterpart of Proposition 5.5 for GUEfe distribution is the following. 
Proposition 5.6. We have 



(5.1) 



E5fc(a;;GUEfe) = exp( ^{xj + 



Proof. Let X be a (fixed) diagonal k xk matrix with eigenvalues xi, . . . ,Xk 
and let A be random k x k Hermitian matrix from GUE ensemble. Let us 
compute 

(5.2) Eexp{Trace{XA)) . 



From one hand, standard integral evaluation shows that (5.2) is equal to the 
right side of (5.1). On the other hand, we can rewrite ( |5.2[ ) as 



(5.3) 



/ PGVE.idy) [ PHaar{du)exp{Trace{YUXU-')) , 

J y\>y2>->yk Ju&u(k) 



where -PcuEfc is probability distribution of GUEfc, Puaar is normalized Haar 
measure on the unitary group U{k) and Y is Hermitian matrix (e.g. diagonal) 
with eigenvalues yi, . . . ,yk. The evaluation of the integral over unitary group 
in (5.3) is well-known, see |Hlj . |H2) . |IZj . |QVj and the answer is precisely 
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Bk{yi, . . . ,yk; xi, . . . , Xk)- Thus, (5.3) transforms into the left side of (5.1 ). 



In what follows we need the following technical proposition. 



Proposition 5.7. Let 



kN 



> 



> 



1,2, 



□ 



be 



a sequence of k- dimensional random variables. Suppose that there exists a 
random variable 0°° such that for every x = (xi, . . . , Xk) in a neighborhood 
of {0, ... , 0) we have 



lim Ei?fe(x; 

A''— !>oo 



kN 



) = ^Bkix; 



Then 



kN 



(j) in the sense of weak convergence of random variables. 



Proof. For k = 1 this is a classical statement, see e.g. |Bi| Section 30]. For 
general k this statement is, perhaps, less known, but it can be proven by 
the same standard techniques as for k = 1. □ 

Next, note that the definition implies the following property for the mo- 
ment generating function of /c-dimensional random variable 4>: 

EBk{xi, . . . ,Xk;a(l) + b) = exp(6(xi H h Xk))^Bk{axi, axu; 4>)- 

Also observe that for any non-constant weakly decreasing f{t) we have 
S{f) > 0. 



Taking into the account Propositions |5.5[ |5.6[ |5.7[ and the observation 
that (e^' — e^i)/{xi 
reduces to the following statement. 



Xj) tends to 1 when Xi,Xj 



0, Theorem 



5.1 



now 



Proposition 5.8. In the assumptions of Theorem 5.1 for any k reals 
hi, . . . ,hk we have: 



lim 

N^oo 



,yNjU) ^ iN-k 



exp 



N^ML{hi + --- + hk) 



Proof. For k = 1 this is precisely the statement of Proposition 4.3 For 
general k we combine Proposition 4.5 and Corollary 3.12 

□ 

5.2. Asymptotics of the six vertex model. Recall that Alternating Sign 
Matrix of size N is a. N x N matrix filled with Os Is and —Is in such a way 
that the sum along every row and column is 1 and, moreover, along each row 
and each column Is and —Is are alternating, possible separated by arbitrary 
amount of Os. Alternating Sign Matrices are in bijection with configurations 
of the six-vertex model with domain-wall boundary conditions. The con- 
figurations of the 6~vertex model are assignments of one of 6 types shown 
in Figure |8] to the vertices N x N grid in such a way that arrows along 
each edge joining two adjacent vertices point the same direction; arrows are 
pointing inwards along the vertical boundary and arrows are pointing out- 
wards along the horizontal boundary. In order to get ASM we replace the 
vertex of each type with 0, 1 or (—1), as shown in Figure^ see e.g. |Kul for 
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more details. Figure [4] in the introduction gives one example of ASM and 
corresponding configuration of the 6-vertex model. 



4* 4* + 



a 




a 




b 




b 




-1 



Figure 8. Types of vertices in the six vertex model divided 
by groups and their correspondence to numbers in ASM 



Let Hn denote the set of all Alternating Sign Matrices of size or, equiv- 
alently, all configurations of six-vertex model with domain wall boundary 
condition. Equip "Jm with uniform probability measure and let cjat be a 
random element of Hat. We are going study the asymptotic properties of wat 
as N ^ CO. 

For "i? G Jtv let ai{'d), Ci{'&) denote the number of vertices in hori- 

zontal line i of types a, b and c, respectively. Likewise, let aj{'&), bj{'&) and 
Cj{-d) be the same quantities in vertical line j. Also let a^^{'d), and 
c*-'('i?) be — 1 functions equal to the number of vertices of types a, h and 
c, respectively, at the intersection of vertical line j and horizontal line i. To 
simplify the notations we view Oj, hi and Cj as random variables and omit 
their dependence on {}. 

Theorem 5.9. For any fixed j the random variable f^"^ weakly converges 



to the normal random variable N(0, y^3/8). The same is true for a^^ 
and ttN-j- Moreover, the joint distribution of any collection of such vari- 
ables converges to the distribution of independent normal random variables 
Ar(0,7378). 



In the rest of this section is devoted to the proof of Thorem 5.9 



6 types of vertices in six vertex model are divided into 3 groups, as shown 
in Figurejsj Define a weight depending on the position (i, j) (i is the vertical 
coordinate) of the vertex and its type as follows: 



q Vj -qUi, c: [q - q)UiVj 
ui, . . . ,un are parameters and from now and till the end 



exp(7ri/3) (notice that this implies q + q 



where i^i, . . . , f at 
of the section we set q 
l;q-q-^ =iV3.) 

Let the weight W of a configuration be equal to the product of weights 
of vertices. The partition function of the model can be explicitly evaluated 
in terms of Schur polynomials. 



N 



Proposition 5.10. We have 
i3e:ijv j=i 

where X{N) = {N - 1, N - 1, N - 2, N - 2, . . . ,1,1,0,0) e GT2n- 
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The following proposition is a straightforward corollary of Proposition 

Em 

Proposition 5.11. Fix any n distinct vertical lines ii, . . . ,in and m dis- 
tinct horizontal lines ji, ■ ■ ■ ,jm o,nd any set of complex numbers ui, . . . , Un, 
vi, . . . ,Vm- We have 



(5.4) EjvH 



k=l 



q - qul\ * 



1-1 



n 

\k=l 



' I «A(iV)(l^^) 



(5.5) E^J] 



-1 2 \ '^i,' 

- qvf^ 



q \j 



q -q 



"31 



n 



i ^A(7V)(12^) 



and, more generally 



(5.6) En I n 



, k=l 



qui 



n 

£=1 



-1 2 \ o-ie 



q -q 



nn 

k=li=l 



{q ^ul-qv'j){q ^-q) 
{q-^ul - q){q-^ - qvj) 



jq n-qui)(q -q) 

(g-l - qul){q-^vl - q) 



-l \ Sx(N){ui, . . ■ ,Un,Vl, . . . ,Vm,l 



2N-n~m\ 



where all the above expectations Ejv are taken with respect to the uniform 
measure onl^- 



We want to study N ^ oo limits of observables of Proposition 5.11 



Suppose that n = 1, m = 0. Then we have two parameters ui = u and 
ii = i. Suppose that as — t- oo we have 

(5.7) u = u{N) = exp{y/VN) 



and i remains fixed. Then we can use Proposition 4.3 to understand the 
asymptotics of the righthand side of (5.4). 



As for the left-hand side of (5.4), note that Cj is uniformly bounded, in fact 
d < 2i because of the combinatorics of the model. Therefore, the factors 
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involving Cj in the observable become negligible as — )• oo. Also note that 
CLi + bi + Ci = N, therefore the observable can be rewritten as 



qe 



N 



'q-^e^y/^ -q 



7-1 



7-1 



qe 



with G satisfying the estimate |lnG(y)| < Cy/yN with some constant C 
(independent of all other parameters). 

Now let z be an auxiliary variable and choose y = y{z, N) such that 



(5.8) 



exp(z/ViV) 



qe 



Now the observable (as a function of z) turns into ( - — q-i-q — 



N 



times 



exp(zai / v N) . Therefore, the expectation in (5.4) is identified with the 
exponential moment generating function for / \/ N. 

In order to obtain the asymptotics we should better understand the func- 
tion y{z,N). Rewrite (5.8) as 



^2,/V]v ^ exp(z/VjV)g-i + g ^ 1 + (exp(z/^/iV) - 1)^ 
q-^+q exp(z/^/iV) 1 + {exp{z/VN) - 1) 

Recall that + q = I, therefore 

2y = \^ (\n{l + q-^{exp{z/VN) - 1)) - ln(l + q{exp{z/VN) - 1))) 

= -iq- q-')z - ^^^z^Vn + ^^-^z^Vn + Oiz^/N) 
Note that the last two terms cancel out and we get 

(5.9) y 
Now we compute 



-zi^ + 0{z^/N). 



qe 



N 



exp 



A^ln 1 



-iV32/v^+0(23Ar-3/2) 



1^ 



exp 



Nqz + qiV3z'^/2 - q^z^/2 + o{l] 



exp 



Summing up, the observable of (5.4) is now rewritten as 



(5.10) 



exp 



Nzi^ - zV2 + o(l) 



exp 



Nqz- z^/2 + o{l] 
as 

at - N/2 



Now combining (5.4) with Propositions 4.3, 4.5 (note that parameter in 
these two propositions differs by the factor 2 from that of (5.4)) we conclude 
that (for any complex z) the expectation of (5.10) is asymptotically 



exp 



A^NyE{f) + AS{f)y' + o{l) 



where / is the function ^-i^- Using (5.9) and computing 

i?(/) = l/4, S(/) = 5/48 
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we get 
(5.11) 



exp 



Nzi 



.^/3 



16 



z' + o{l) 



Now we are ready to prove Theorem 5.9 



Proof of Theorem 5.9 Choose Zk and z'^ to be related to Uj, and Vj., re- 



spectively, in the same way as z was related to u (through (5.7) and ( |5.8| )). 
Then, combining the asymptotics (5.11) with Corollary 3.12 we conclude 
that the righthand side of (5.6) as — )• oo is 
(5.12) 



JJexp 



fc=i 



Nzki 



16 



4 + 0(1) 



JJexp 



Nz'd 



16 



(4r + o(i) 



Now it is convenient to choose Zi (z^) to be purely imaginary Zi 

Summing up the above discussion, observing that the case n = 0, m = 1 
is almost the same as n = 1, m = 0, (only the sign of aj changes), and that 
the observable (5.6) has a multiplicative structure, and the third (double) 
product in (5.6) is negligible as N ^ oo, we conclude that as iV — )• cxd for 
all real Sj, s' 



(5.13) lim Eat exp 

N-^oo 



^au-N/2 . ^ai,-N/2,. 



.k=l 



N 



exp 



3 

16 



\k=l 



The remainder o(l) in the left side of (5.13) is uniform in aj^., and, 
therefore, it can be omitted. Indeed, this follows from 



Etv exp 



ai-N/2 ■ 
-SI + o(l) 



■ E TV exp 
<Etv 



ai - N/2 ; 

SI 



exp 



ai - N/2 ; 

^ SI 



o(l) = o(l). 



Hence, (5.13) yields that the characteristic function of the random vector 



ai^ - N/2 ai^ - N/2 aj, - N/2 



N 

converges as — )• oo to 
exp 



A^ 



N/2 \ 



N J 



3_ 
16 



K,k=l e=i y 

Since convergence of characteristic functions implies weak convergence of 
distributions (see e.g. [Bil Section 26]) the proof of Theorem 5.9 is finished. 

□ 
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5.3. Towards dense loop model. In |GNP] de Gier, Nienhuis and Pon- 
saing study the completely packed 0(n = 1) dense loop model and introduce 
the following quantities related to the symplectic characters. 
Following the notation from |GNP] we set 

tl(zi, ...,zl) = X\l{zI, ...,zl) 

where A-^ e GT^ is given by Xf = [^^J for i = 1, . . . , L. Further, set 
(5.14) 

TL+l(Cl,^l, • • ■,Zl)tl+i{C2,Zi, . . .,Zl)' 



^^l(Ci, C2; zi, . . . , zl) = (-l)^i— In 
Define 



Tl{zi, Zl)tl+2{Ci,C2, Zl,..., Zl) 



and 



(7) 

De Gier, Nienhuis and Ponsaing showed that Xf^ and Yl are related 



= ZjT—UL{Cl,C2;Zi,...,ZL) 
OZj 

d 

Yl = W — UL+2{Cl, C2; Zl,..., ZL, vq~^,w)\y,=v, 
(i) 

in particular, Xf^ is a function oi zi, . . . , zl and Ci) 5 while Yl also depends 
on additional parameters v and q 
De Gier, Nienhuis and Ponsai 
to the mean total current in the 0{n = 1) dense loop model, which was 
presented in Section 1.6 More precisely, they prove that under certain fac- 
torization assumption and with an appropriate choice of weights of configu- 

ii) 

rations of the model, Xf^ is the mean total current between two horizontally 
adjacent points in the strip of width L: 

and Y is the mean total current between two vertically adjacent points in 
the strip of width L: 

Yl = F^^''^^'^^~^''^\ 

see [(xNPj for the details. 

This connection motivated the question of the limit behavior of Xj^ and 

as the width L tends to infinity, this was asked in |Gi| . |GP| . In the 
present paper we compute the asymptotic behavior of these two quantities 
in the homogeneous case, i.e. when Zi = 1, i = 1, . . . , L. 

Theorem 5.12. As L ^ 00 we have 



and 



Y, 



=1, j=i,...,L 4L \L 



(7) 

Remark 1. When z = 1, Xf^ is identical zero and so is our asymptotics. 
Remark 2. The fully homogeneous case corresponds to w = exp"''^/^, 
q = e^'^'/'^. In this case 
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Remark 3. The leading asymptotics terms do not depend on the bound- 
ary parameters Ci and ^2- 

The rest of this section is devoted to the proof of Theorem 5.12[ 

Proposition 5.13. The normalized symplectic character for = 
( , 5 • • • > li 0, 0) is asymptotically given for even L by 



Se-h^ey-l) f i (e3/2s/-l)2 



{e^/^y - i){ey + 1) \ 9 ey/^ey - 1)2 

and for odd L by 



3e--4y{ey - 1) (a {e^i'^y - 1) 



(e3/2?/ - \){ey + 1) I 9 evl'^{ey - 1)2 
for some analytic functions ti, t2, ... and t'^^t^, ■ 



1 I h{y) . h{y) 

Li/2 ^ L2/2 ^ • ■ ■ 



LV2 ^ L2/2 ^ ■ ■ ■ 

such that ti = t'l and 



t'2 = t2 + ^^{e'/'y 



lYe 



-3/2y 



Proof. We will apply the formula from Proposition 3.19 to express the nor- 
malized symplectic character as a normalized Schur function. The corre- 
sponding v is given by Ui = + 1 for i = 1, . . . , L and Ui = — 
for i = L + 1,...,2L, which is equivalent to = + 1 for all 
i = 1, . . . , 2L. We will apply Proposition 4.2 to directly derive the asymp- 
totics for Sy{ey\ 2L, 1). For the specific signature we find that f{t) = | — 5* 
and 

Hw;f) 



[ ln{w - f{t) - 1 + t)dt 
Jo 



_ 1 
~ 6 

In particular, we have 



-6 + {5-4w)ln 



+ w 



+ (1 + 4w)ln 



+ w 





5 


-In 


"1 


-- +w 


- +w 




4 




4 



:F"iw;f) 



1 



The root of J-'{w; f) = y, referred to as the critical point, is given by 

1 + 5e3/22/ 



wo = wo{y) 



4(-l + e3/23/)' 



Example 2 of Section 4.2 shows that a steepest descent contour exists for 
any complex values of y for which wq [—1/4,5/4], i.e. if e^/^^ is not a 
negative real number. The values at wq are 



and 



ywo - F{wo-f) = - 
J"'{wo;f) 



+ ln{e^/^y - 1) + 1 - In ^ 



4 (e3/2i/ _ 1)2 

'9 e^/^y 
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In order to apply Proposition 4.2 we need to ensure the convergence of 



as 



Q{w] u, /), defined in Section 4.1 

ni n (, I /{j/{2L))-,^,/{2L) \ 

Qiv,. /) = E (1 + .-/0/(2L))-l+,/(2Z,)j ■ 

Substituting the values for using the formula ln(l + ~ x, and approxi- 
mating the sums by integrals we get 

(5.15) Qiw- /) = E ^ _ _ 1 + ,/(2L) + = 

2L (L-i)/2+l 1 I i 2L (L-i)/2+l/2 1 1 » r^ 

E 2L "'" 4 2 2L I 2L "'" 4 2 2L , ^V"^; ^) 

i=l,feL(mod2) 4 + 4L ^ 2L i=l, i=L+l(mod 2) 4 + 4L ^^^21 

2L _J_ 2L _J_ Af r\ 

i=l,feL(mod 2) 4 ^ 4L j=l, i=L+l(mod 2) 4 ^ 4L 



5 , 3^-^'?+ / ^n57^^ + + r + 



ou;-| + fr/ Jo w-l + ^r] L L 

-2^"l^^. + lJ+ L + L • 

for some functions yl(t(;;L), B{w;L) bounded in w and L. In fact we need 
to understand how A{w;L) and B{w;L) depend on the parity of L. First, 
note that (using ln(l + x) w x — x^/2) 



The last sum can be again approximated by the integral similarly to (5.15); 
therefore 

A{w;L) = A{w) + 0{l/L). 

Next, B{w, L) appears when we approximate the integral by Rieman sum. 
Since trapezoid formula for the integral gives 0(1/L^) approximation, we 
have (with the notation u{x) = ——. — \ — 5-^) for even L 

B{w; L) = -n(0) + u + 0{l/L) = u{l) - m(0) + 0{l/L) 

and for odd L 

B{w; L) = -u(0)/2 + u /2 + 0{l/L) = u{l)/2 - u{Q)/2 + 0{l/L) 
Therefore, we have 

A{w, L) + B{w, L) = C{w) + (-1)^+^^ {-^ - + 0(1/ L) 



ASYMPTOTICS OF SYMMETRIC POLYNOMIALS 



53 



and hence we obtain as L — )• oo 



exjp{Q{w;u,f)) 
and 



w — 
w + 



5 \ 2 

4 



i + (-i) 



L+l 



16L 



W - i U7 + 4 



exp(Q(^o;^,/)) = exp(-^y) (^1 + (-1)^+'^ {{e'^'' " l)'e-=^/2^) + 0(1/l2)) 

Substituting into Proposition 4.2 the expansion of Q and exphcit values 
found above we obtain 



(5.16) S^{ey;2L,l) 



wo - /(O) - 1 



T"{wo){wo-f{l)) Uo + i 



exp 2L{ywQ - F{wq)) 



X l + (-l) 



L+l 



1 



1 



16L 



+ ... (1 + ...) 

^«o - i w^o + 4 / / 



3e-!s/(e!^ - 1) A (e3/2y - 1)2 



2(e3/2?/_i) I 9ey/2(e?/ - 1)2 



1 + U-1/2+ t2 + (-l)^+' 



(e3/2y _i)2g-3/2y^ 



24 



Proposition 3.19 then immediately gives X_)^L(e^; L, 1) as -^^Sy{ey]2L,l). 

□ 

We will now proceed to derive the multivariate formulas needed to com- 
pute ul- First of all, notice that 3i\{x] L) = aL(x)/i(a;)^^^(2 — x — x~^)~^ 
for h{x) = |3;-3/2(x3/2_i)2 and aiix) = a{x)+bi{x)L~^/'^ +b2{x)L-^ + . . . . 
The final answer does not depend on the precise formulas for a and bi . Define 

~ r \ X\L{zf,...,Z^,l^ ^) ^ / 2 2 r i\ 

TL{zi,...,Zk) = T-FT =X;^L(Zi,...,Zfc;L,l) 



ul{C,i,C,2]zi,. . . ,zk) = (-l)^i 



V3 



In 



TL+l{Cl,Zl, . . . ,Zk)TL+liC2,Zi, ...,Zk) 



TL^Zl, . . . ,Zk)TL+2iCl,C2,Zl, . . . , Zk) 

Then ul{Ci, C2] zi, . . . , Zk) — ul{Cii C2, zi, . . . , Zk) is a constant and thus we 
have 

d ~ ^(i) 



^i^S^(Cl,C2;^l) = x^^^ 



d 



w-^UL+2{CiX2;vq ^,w) 



Yl. 



Therefore, we can work with X\l instead of x\l and with u instead of u. 
For any function ^ and variables t>i , . . . , Vm we define 
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Proposition 5.14. Suppose that 

X\{x;L, 1) = ql(x)/i(x) 

where 



aiix) = a{x){l + bi{x)L ^^'^ + b2{x)L ^ + ...) for even L, 

aL{x) = a{x){l + 6i(x)L"^/2 + h2{x)L-^ + . . . ) for odd L 

and a{x),bi{x),b2{x), h{x) are some analytic functions of x and let ^(x) 
x-^ ln{h{x)). Then for any k we have 

Xa(xo, ...,Xk;L + 1)Xa(xi, . . . , Xk+i, L + 1) 



In 



Xa(xi, . . . , Xk; L)Xx{xo, Xk+r, L + 2) 



ci(xo,Xfe+i;L) + ^2(62(2;) - h2{x)) — = h 



i=l 



L 



In {i{xk+if -iixof) + -(B(xo,..,Xfc;0 - B{xi, ..,Xk+i;Cj + C2(xo, x^+i)) 

Proof. Apply Theorem |3.17| to express the multivariate normalized character 
in terms of OLL{xi) and /i(xj) as follows 



(5.17) 



2ix{xi,. . . ,Xm;N) 



\{^X{X^■,N) 

i=o 



(2iV- 1)! 



(^1 ; • • • ) 2^71 



which is applied with A^ = L,L + l,L + 2, m = A;,/c + l,/c + 2 and define for 
any N and m 



(5.18) M7v(xi,...,x^) :=det 



Z^f^'-^ K(xO/i(x,)^] 
N^i-^aN{xi)h{x^)^ 



ni^i aAr(2;j)/i(xi)^ 



ni!Li aN{xi)h{xi 



where, as above, Di = XiJ^. The second form in (5.18) will be useful later. 



We can then rewrite the expression of interest as 

"Xa(xo, . . . , Xfc; L + 1)Xa(xi, . . . , x^+i; L + 1) 



(5.19) In 



+ In 



n 



Xx{xi, . . . , Xfc; L)Xa(xo, . . . , Xk+i;L + 2) 

Xa(xo; L + l)XA(xfc4.i; L + 1) 
Xa(xo; L + 2)XA(xfc+i; L + 2) 

In 



= consti{L) + In 



XA(xi;L + l)^ 



L Xx{xi;L)Xx{xi;L + 2) 



(Xo - 1)^X0 ^(Xfc+i - ifx ^ 



k+l 



Xq + Xq^ - (Xfc+i + X^,^^) 



In 



Ml+i(xo,xi, . . . ,Xfc)ML+i(xi, . . . ,Xfc+i) 

ML{xi,...,Xk)ML^2{xQ,...,Xk+i) 
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where consti{L) will be part of ci{xo, Xk+i] L). We investigate each of the 
other terms separately. First, we have that 



Xx{x;L)Xx{x;L + 2) 



E 

i=l 



In 



aL+ijxi)'^ 
aL{xi)aL+2{xi] 



+ ln 



aL+i{xo)aL+i{xk+i 
aL+2{xo)aL+2{xk+i) 



where the terms involving xq and x^+i are absorbed in ci and we notice 
that 



In 



aL+iix)' 



Next we observe that for any i and 



2(b2{x)-b2ix))- 



-1 



L 



+ 



1 

Z2 



(5.20) 



[x-^] [aN{x)h{x 
N^aN{x)h{x)N 



i{xf + 



qi 



CixY + iri^ix) 



log(a(x)). In particular, 



where qi = ^ix){x-^£,{x) + ^{x^) and ri{x) = x-^ 

since M^r is a polynomial in the left-hand side of (5.20), it is of the form 
(5.21) 

Mn{xi, ...,Xm)= Mfixi), . . . , fiXrr.)) + Pl (xi , . . . , X„) - + O 

for some function pi which depends only on ^ and a. That is, the second 
order asymptotics of Mj\f does not depend on the second order asymptotics 
of ol- Further, we have 



Mn 

Mn+] 



1 + O ( N-^/"^ 



for any N, so in formula (5.19) we can replace M^^+i and Mi_|_2 by Ml 
without affecting the second order asymptotics. Evaluation of M directly 
will not lead to an easily analyzable formula, therefore we will do some 
simplifications and approximations beforehand. 

We will use Lewis Carroll's identity (Dodgson condensation), which states 
that for any square matrix A we have 

(det A)(det Ai,2;i,2) = (det Ai;i)(det^2;2) - (det Ai;2)(det ^2;i), 

where Aj-^j denotes the submatrix of A obtained by removing the rows whose 
indices are in / and columns whose indices are in J. Applying this identity 
to the matrix 



A 



L'^^aL{xi)h{xi)^ 



k+l 



ij=0 



we obtain 



(5.22) ML{xi,...,Xk)MLixo 



det 



Dl'{aL{xi)h'^{xi)) 
L'^^aL{xi)h{xi)^ 



i = [0:fc-l,fc+l] 



det 



[l:fc+l] 



D^'{aL(xi)h^{x.i)) 
L'^^aL{xi)h{xi)^ 



i,j=0 
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det 



L^iaL{xi)h{xi)^ 



j=[0:k-l,k+l] 



det 



i=[0:k] 



Dl\aL{xi)h^{xi)) 
L'^iaL{xi)h{xi)^ 



j,i+l=l 



where [0 : A: — 1, /c + 1] = {0, 1, . . . , A; — 1, A; + 1}. The second factors in 
the two products on the right-hand side above are just Ml evaluated at the 
corresponding sets of variables. For the first factors, applying the alternate 



formula for Ml from (5.18) and using the fact that 



/\{vi, ...,Vm)^^Vi = det 

i=l 



j = [0:m— 2,m] 
i=[l:ra\ 



we obtain 



det 



Dl^{aL{xi)h^{xi)) 
L'^iaL{xi)h{xi)^ 



j = [0:fc-l,fc+l] 



j=[l:fe+l] 



det [(AVi^^)^]dl:ti''^'^ n O^L{x.)Kxd' 

i = l 

k-\-\ \ 



k+l 



Ui=lC^L{x^)hix,)L^^" Lv--./^ ; U=[V.k+l] 



Ui=i aL{xi)h{xi)^ \ 

1 /fc+i 



Ki=l 



^k+1 



vj = l 



WtZ{<^L{Xi)h{Xi, 

Substituting these computations into ( 5.22[ ) we get 

Ml(xo,Xi, . . .,Xk)ML{xi, . . .,Xk,Xk+l) 



Jl aL{xi)h{xi)^ Ml{xi, . . .,Xk+i) 



(5.23) 



Ml{xi, . . .,Xk)MLixo 



?j) (nr^i «L(x.)/i(x.)^)Mi(xi, . . . , x,+i] 



rfc+l 



lli^l aL{xi)h{xi)^ML{xi, Xk+i) 



(Ej=o Tj) (ULo aL{xi)h{xi)^)ML{xo, ...,Xk] 



nj=o aL{xi)h{xi)^ML{xo, ■..,Xk) 
Dl^^aLixk+i)h{xk+i)^ DlaL{xQ)h{xQ) 



+ 



L'^aL{xk+i)h{xk+i)^ L'^aLixo)h{xo)^ 



L'^Ml{xi, . . .,Xk+i) 

Di[aLixi)h{xi)^] DiMLjxi, x^+i) 
^ LaL{xi)h{xi)^ LMl{xi, . . . ,Xk+i) 



L'^Ml{xq, ...,Xk) 



E 

i=0 



Di[aLixi)h{xi)^] DiMLjxo, ■ ■ . ,Xfc) 
LaL{xi)h{xi)^ LMl{xq,. . . ,Xk) 



Using the expansion for Ml from equation (5.21) and the expansion from 



(5.20) we see that the only terms contributing to the first two orders of 



approximation in (5.23) above are 



1 



(5.24) i{xk+lY - i{xoY + -(C3(xfc+i) - C3(xo)) + 



ASYMPTOTICS OF SYMMETRIC POLYNOMIALS 



57 



2 
L 



^^^^^^ Aiax,r,...,axk+,r) ^^^""^ 



j=0 



A(e(xo)2,...,^(:Efc)2) 



for some function C3 not depending on L, so C2(xo,Xfc+i) = C3(xfe+i) — C3(xo). 
Substituting this result into (5.19) we arrive at the desired formula. 

□ 



Proof of Theorem 5.12 . Proposition 5.14 with xq = Cl; ^k+i = Ci and Xi 
zf shows that 

(5.25) l[ul{(:iX2,zi,...,zu)- ci(C?, C|; L) - ^ 2(62 (x) - 62 (x)) ^ 



i=l 



-In 



(C(C?)' - e(C2)') + 2(i?(xo, .., Xfc; e) - .., Xfc+i; e) + C2(Ci', C2)) 



>2\2 



2 /'2\ 



converges uniformly to and so its derivatives also converge to 0. Proposi- 
tion [STlS] shows that in our case 

h{x) = ^X-3/2(x3/2 _ 1)2 



and thus ^(x) 
equation 



9 

i ■ fsTT^- Moreover, the function ^ satisfies the following 



9 x^/^ 



9 



(e(x)2 - 1) 



2(x3/2-l)2 8 

and so we can simplify the function i? as a sum as follows 



(5.26) Bivi,...,vm;0 



E.^(^0^»^A(g(t;i)2,...,^(^„ 
A{avir,...,avm)^) 



Vi Vi 



EE 



Z^Z^ ^(^.)2_^(„)2 



e(^,)2-C(t;,)2 



= E-i(^(-^)'+^( 

i<j 

We thus have that 



1) 



(m-1) (^^(^^0') + 



9 j m 
412 



S(xo, . . . , Xfc; ^) - 5(xi, . . . , Xfc+i; ^) 
which does not depend on xi, . . . , x^,. 



9 



mxk+if-a^of), 



Differentiating (5.25) we obtain the asymptotics of Xf^ as 



X 



.V3 



d 



{-l)^Z — 2(b2{z^)-b2{z'))- 



-1)^ .\/3 d 



L 



1 z- 

2 dz 



{z-' - Ifz 



2^-3 



.V3 3 _q> 

) 
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(7) 

For Y£ the computations is the same. 



□ 



6. Representation-theoretic applications 

6.1. Approximation of characters of U{oo). In this section we give a 
new proof of Theorem 1.5 presented in the Introduction. 

Recall that a character of U{oo) is given by the function x{ui,U2, ■ ■ ■), 
which is defined on sequences Ui such that Ui = 1 for all large enough i. 
Also x(l)l5---) = 1- By Theorem 1.3 extreme characters of C/(oo) are 
parameterized by the points uj of the infinite-dimensional domain 

n C M^~+2 = X M°° X M°° X M°° X M X M, 

where O is the set of sextuples 

UJ = (a+,a-,/3+,/3-;(5+,(5-) 

such that 

a± = (af > > • • • > 0) G M°°, = > > • • • > 0) G M°°, 
00 

i=l 

Let /i be a Young diagram with the length of main diagonal d. Recall 
that modified Frobenius coordinates are defined via 

Pi = Hi-i + 1/2, Qi = fi'i - i + 1/2, i = l,...,d. 

Note that Yli=iPi + Qi = l/^l- 

Now let A G GTtv be a signature, we associate two Young diagrams A+ 
and A~ to it: row lengths of A"*" are positive of Aj's, while row lengths of A~ 
are minus negative ones. In this way we get two sets of modified Frobenius 
coordinates: pf ,qf , i = 1, . . . ,d'^ and p~ ,q~ , i = 1, . . . ,d~ . 

Proposition 6.1. Suppose that X{N) G GTjv is such a way that 



^^a+ ^^/3+ ^^a- ^ 



i=l i=l 



then 



TV- 

where 



lim S'wAr)(x; A, 1) = $00 ( a,/3,7; 



$00 (^a,^,7;;^^) =exp(7+(x-l) + 7 (x ^ - 1)) 

l + /3+(x-l) l + (l-/3-)(x-l) 



00 



T/ie convergence is uniform over 1 — e < \x\ < 1 + e for certain e > 0. 
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Remark 1. Note that 

l + (l + a-)(x-l) l-Qi(x-i-l)^ 



which brings the functon <I>oo into a more traditional form of Theorems 1.3 



Remark 2. Our methods, in principle, allow us to give a full description 
of the set on which the convergence holds. 

Proof. The following combinatorial identity is known (see e.g. [BO, (5.15)] 
and references therein) 

(61) yr S + i-K _yr s + 1/2 -pf + 1/2 + N+pf 
t\ « + ^ l\s + l/2 + qtl\s + l/2 + N-qr 

Introduce the following notation: 

d+ i/^-pf d- 1/2+N+P+ 

^ " ' ii l/2+g+ 11 , l/2+N~q- 



i- 



--IW + .J ' i=l W + 



and observe that ( |6.1[ ) implies that (here f{t) = 0) 

(6.2) exp {Q{w; X{N), /)) = $^(A(iV); u;). 
As — 7- oo we have 

(6.3) ^N{X{Ny,w) ^ «>oo(a, /3, 7; w) 

uniformly over 1— ffi < \w\ < 1 + ei for some ei > 0. In fact we have the 
convergence on much greater set, which is, basically, all C except for the 
poles of $00 • 

Now we can use Propositions 4.2 and 4.4 with the steepest descent con- 
tours of Example 1 of Section |4.2[ Recall that here f{t) = 0, x = e^, 
J'{w;0) = wln{w) - (w - l)ln{w - 1) - 1 and wq = 1/(1 - e"^). We 
conclude that as A" — )• 00 



5aw (e^; N,l) = ^ -r'{wo;0){wo) ^^^^^ e^iev - 1)^-^ ['^'^'^ ) ' 

Substituting wq, Q{w; X{N), f), using (6.2), (6.3) and simplifying we 
arrive at 

Sx(Ny,N,i{^) (^a,/3,7; ^^^^ 

The above convergence is uniform over compact subsets of 

V = {x = e C \ -TT < Im{y) < vr, £2 < Re{y) < £2 2/ / 0}. 

(Here the parameter £2 shrinks to zero as af goes to infinity.) 

It remains to prove that this implies uniform convergence over 1 — e < 
|x| < 1 + e. 

Decompose 



x\ 
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Since Sx(n) is a polynomial, only finitely many coefficients Ck{N) are non- 
zero. The coefficients Ck{N) are non-negative, see e.g. [Ml Chapter I, Section 
5], also Zk Ck{N) = SxiN){l; N, 1) = 1. 

Since <I>oo(a, /3, 7; ;jir[) is analytic in the neighboorhood of the unit circle 
we can similarly decompose 



a, (3 



1 \ °° 



(6.4) 



CkiN) 



Sx(N)(.z;N,l)z 



-k-l 



dz 



k=--oo 

We claim that liuii^^oo Ck{N) = Cfc(oo). Indeed this follows from the 
integral representations 

and similarly for <I>oo- Pointwise convergence for all but finitely many points 
of the unit circle and the fact that \Sx{n){^j ^ 1 for \z\ = 1 implies 

that we can send — )• 00 in (6.4). 

Now take two positive real numbers 1 — ei<a<l<6<l-|-ei such that 

(6.5) lim SxiN) («; N, 1) = (a, /3, 7; ) , 



(6.6) 



lim Sx^N)ib;N,l) 



a,/3,7; 



1 



1 



For w satisfying a < \w\ <b and some positive integer M write 
(6.7) 



5^(cfc(iV) - Ck{oo))w' 



< j;|cfe(iV)-Cfc(oo)|(a'= + 6^ 



AI 



< 



|cfc(iV)-Cfc(oo)|(a'=+6'=)+ Y Ck{N){a'^+b')+ ^ Cfe(oo)(a'=+6'=) 

fc=-Af |fc|>Af |A:|>M 

The third term goes zero as M — )• 00 because the series ca;(oo)z'^ con- 
verges for z = a and z = b. The second term goes to zero as M — )• 00 
because of (6. 5), (6. 6) and Ck{N) — Cfc(oo). Now for any 5 we can choose 
M such that each of the last two terms in ( |6.7| ) are less than 6/3. Since 
Ck{N) — )• Cfc(oo), the first term is a less than 6/3 for large enough A^. There- 
fore, all the expression (6.7) is less than 6 and the proof is finished. □ 



Now applying Corollary 3.10 we arrive at the following theorem. 



Theorem 6.2 (cf. Theorem 1.5). In the settings of Proposition 6.1 for any 
k we have 

lim Sx(N){xi, . . . ,Xfc;iV, 1) = TT^oo (a,P,T, — ^ ) , 

N->oo fJl \ Xi-lJ 

The convergence is uniform over the set 1 — £ < \xi\ < 1 + £, i = 1, . . . ,k 
for certain e > 0. 
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Note that we can prove analogues of Theorem |1.5| for infinite-dimensional 
symplectic group Sp{oo) and orthogonal group 0(oo) in exactly the same 
way as for U{oo). Even the computations remain almost the same. This 
should be compared to the analogy between the argument based on binomial 
formulas of [OOj for characters of U{oo) (and their Jack-deformation) and 
that of |002j for characters corresponding to other root series. 

6.2. Approximation of (7 deformed characters of [/(oo). In [GJ a q- 
deformation for the characters of [/(oo) related to the notion of quantum 
trace for quantum groups was proposed. One point of view on this deforma- 



tion is that we define characters of [/(oo) through Theorem 1.5, i.e. as all 
possible limits of functions 5a, and then deform the function Sx(n) keeping 
the rest of the formulation the same. A "good" ^-deformation of turns out 
to be (see [G] for the details) 

sxjxi, ...,Xk,q ,-..,q ) 

Throughout this section we assume that q is a real number satysfying < 
q < 1. The next proposition should be viewed as g-analogue of Proposition 
EH 

Proposition 6.3. Suppose that X{N) is such that A7v_j_|_i — t- for every 
j . Then 

sx{x,q~'^,q~'^, . . .,q^~^) ^ j-, . . 

F,y[X), 



sx{l,q-\...,q^-^) 
L - qi+^) Injq 



(6-8) F^i^) = ft / ——^ —dz, 



j=0 

where the contour of integration C consists of two infinite segments of 
Im{z) = 9'^'''^9 right and vertical segment [— M(C') — 

jj^,— M(C') -|- i^^] with arbitrary M[C') < ui. Convergence is uniform 
over X belonging to compact subsets o/C \ {0}. 

Remark. Note that we can evaluate the integral in the definition of 
Fy{x) as the sum of the residues: 

(6.9) Fu{x) = W\ ^—lYl-, 1 yT\ ^-^Kdz, 



The sum in (6.9) is convergent for any x. Indeed, the product over j > k 
can be bounded from above by l/(q;q)oo- The product over over j < k is 
(up to the factor bounded by {q]q)oo) 

k-l 

JJ q^k+k-yj-j ^ 
i=l 

Note that for any fixed m, if /c > ko^m) then the last product is less than 

qm{uk+k-l) 



62 



VADIM GORIN AND GRETA PANOVA 



We conclude that the absolute value of A:th term in (6.9) is bounded by 

I |!/fc+fc-l m(i/fc+fc-l)_ 



1 



Choosing large enough m and k > ko{m) we conclude that (6.9) converges. 



Proof of Proposition 6.3. We start from the formula of Theorem 3.6 
(6.10) 



s\{x, l,q~ 



J2-N\ 



1-N\ 



n 

i=0 



l-N 



Hq) 



{x- 



2tti 



{x/qY 



cuU 



-dz 



where the contour contains only the real poles z 
the rectangle through M + M - -M - 
sufficiently large M. 
Since 

sx{x,l,q~^,. 
we may also write 
(6.11) 



Xj + N - j, e.g. C is 
-M + 4) for a 



Tri 
ln(q) 



)=qWsx{q-'x,, 



l-N 



sx{x,q ^,q ^) 



^-2 t^l-N 



l-N] 



n 



(gi 



Hq) 



{x 



N-2 



(1 



j=0 

-g^+^)ln(g) 



27ri 



j=0 



2tt\ 



cnf=i(9-^-9 



-A, 



c^n =1(1-9- v^+^-^) 



Note that for large enough (compared to x) N the integrand rapidly 
decays as Re{z) — >• +00. Therefore, we can deform the contour of integration 
to be C which consists of two infinite segments of Im(z) 



-H(g) going to 

the right and vertical segment [— M(C') — j^^, —M{C') + jjf^] with some 
M(C'). 

Note that prefactor in (6.11) converges as — ?■ 00. 

Let us study the convergence of the integral. Clearly, the integrand con- 
verges to the same integrand in F^{x), thus, it remains only to check the 
contribution of infinite parts of contours. But note that for z = s i: 
s G M, we have 



Hi) 



X 



nf=i(i-«-^g^^-^^-^-) nf 



:i+^ 



jX,+N-j-^ 



Now the absolute value of each factor in denominator is greater than 1 and 
each factor rapidly grows to infinity as s — t- 00. We conclude that the 
integrand in (6.11) rapidly and uniformly in N decays as s — >• +00. 



It remains to deal with the singularities of the prefactors in (6.11) and 



(6.12) at X = q~^. But note that pre-limit function is analytic in x (indeed 
it is a polynomial) and for the analytic functions uniform convergence on a 
contour implies the convergence everywhere inside. □ 



As a side-effect we have proved the following analytic statement 
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Corollary 6.4. The integral in (6.12) and the sum in (6.9) vanish at x = 
q-\ 

Uj for every j. 



Theorem 6.5. Suppose that X{N) is such that X^^j^i 
Then 

- ^FW(xi,...,Xfc), 



•5A(A)(1 



(6.12) 



(-l)©g-2(3) 



k 



det [D^r^l,] l^^^ n Fu{x^q''-^){xq^-^■ qU 
1=1 

Convergence is uniform over each Xi belonging to compact subsets o/C\{0}. 



Remark. The formula (6.12) should be viewed as a g-analogue of the 
multiplicativity in the Voiculescu-Edrei theorem on characters of U{oo) 
(Theorem 1.3). There exist a natural linear transformation, which restores 
the multiplicitivity for g-characters, see [G] for the details. 



Proof of Theorem\6.5\ Using Proposition |6.3| and Theorem 3.5 we get 



F«(xi, ...,Xk)= lim g-'^l^(^)ISA((z'xi, . . . , q'^x^; N, q-') 

A— !>oo 

{k+i\ 



X 



(-1)© det [D=-%^, ^ S^{x,q^- N, q-^) nfj^x.g'^ - r^+^) 



[N-iL-i\ 

In order to simplify this expression we observe that 7 ^ r — )• 

[iV-l],-i! 

m 

^Nii-i)-{"/) as AT ^00. Also, Ylixq'' - q-^+^) = {-l)'^q-^^^) {xq^-^]q)m- 

i=i 

Last, we have 

lim q~\^\Sx{xq^]N,q-^) = F^{q''-^x). 
A-s>oo 

Substituting all of these into the formula above, we obtain 
FW(xi,...,Xfc)= hm ^- fci^ X 



q 



1 (-l)®det[l?-l,]^^^,^ 



llF,{x,q'^-'){xq'^-';q)c 
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